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INTRODUCTION 


Purpos 


The  Columbia  River,  which  drain:'  a  large  portion  of  the  Pacific 
northwest,  discharges  a  mean  annual  av<  rage  f  .  57,000  ft“  sec  '  (7,300  11 

)  f  fresh  water  into  the  northeast  Pacific  Ocean  at  U6°  20'N.  This 
amounts  tc  lU  percent  of  the  total  annual  discharge  from  the  United  States 
f  Ameri  'a  as  estimated  by  Langbein  (19^9)*  The  flow  is  seasonal  in  nature 
with  a  maximum  in  June  and  a  secondary  maximum  in  January.  I:  never  dr  t . 
bel  w  approximately  half  of  its  average  rate  and  the  River  at  all  seasons 
•ontributes  the  majority  of  fr  shwator  discharged  from  the  Pacific  coast 
between  southern  California  and  southeastern  Alaska. 

The  area  of  the  Pacific  Ocean  into  which  the  ColumDia  River  dis¬ 
charges  is  -hara  :terizcd  by  weak  and  poorly  defined  currents  (Barnes  and 
Paquette  1957).  The  west-wind  drift  of  the  North  Pacific  Ocean  sets  east¬ 
ward  and  on  approaching  the  continent  of  North  America,  diverges  well 
offshore  near  U5JN.  The  northern  oranch  feeds  a  gyre  in  the  Gulf  of  Alaska 
and  : he  southern  branch  turns  southward,  forming  the  California  Current. 

The  Columbia  -ntero  the  Pacific  from  a  straight  North-South  coast  line  and 
•an  disperse  through  an  arc  of  180  degrees. 

The  low  salinity  water  associated  with  the  river  discharge  is  dis¬ 
tributed  in  the  form  of  a  plume.  This  plume  extends  from  the  river  mouth 
over  many  tens  of  thousands  of  square  miles  and  remains  identifiable  for 
long  periods  of  time.  The  discharge  of  the  Columbia  River  into  the  north¬ 
east  Pacifi  •  therefore  provides  a  good  system  for  the  study  of  the  nature, 
movement,  dispersion  and  related  effects  of  river  effluent  in  the  ocean, 
as  well  as  a  very  good  natural  system  for  studying  the  effects  of  wind  on 
ocean  surface  waters. 

J 

The  Department  of  Oceanography,  University  of  Washington,  under  con¬ 
tract  with  the  Atomic  Energy  Commission  from  1961  to  the  present,  has 
intensified  its  studies  in  the  area  of  Columbia  River  discharge.  The  first 
objectives  were  to  determine  the  gross  features  of  the  movement  and  dis¬ 
persion  of  the  effluent  water  in  the  open  sea,  and  ultimately  to  ascertain 
the  probable  fate  of  any  radioactive  material  that  might  be  associated  with 
river  effluent  or  river-borne  materials.  Studies  of  the  physical,  chemical 
geological  and  biological  properties  of  the  effluent  water  and  ambient 
ocean  water  were  undertaken  with  these  objectives  in  mind.  This  paper 
reports  the  gross  features  of  the  plume  as  they  were  observed  in  1961  and 
19o2  and  describes  quantitatively  certain  features  of  the  dispersion  of  the 
plume  water;  other  results  of  this  continuing  study  will  be  reported  when 
analyses  are  completed.  The  intense  sampling  program  which  commenced  in 
January  1961  has  formed  the  basis  for  investigations  of  chemical  properties 
by  Richards  and  Stefansson  (1963)  and  Stefansson  and  Richards  (196^); 
biological  studies  by  Anderson  (196^);  and  geological  studies  by  Gross  and 
co-workers  (1963). 


Th'*  observation;  and  descriptions  presented  in  this  pap  -r  arc  based 
n  the  1  *  .  and  .  .  r  sulti  f  f  truises  of  the  R.  V.  Brown  Bj  ir  1 

area  f fluent  f  th  ’  I  La  River,  Pwe Lvf  f  thes< 

cruise;  w -re  designed  specifically  to  delineate  tin*  f  effluent  dis¬ 

perse  n.  Jrui  ■  sine  1962  ha\  investigated  urn  1  th<  fin  r  ietail. 

I:;  addition  to  th-  n.-n; areni'-nts  of  temperature,  salinity,  and  exjv.  n,  made 
r  ttinely  al  ach  :eai  raphic  station,  :ontinuou!  measurement  >f  trface 

salinity  and  temperature  by  means  of  a  salinity-t emperutur--  deviation 
••  ■  rder  were  made  while  underway.  At  many  stations  .  tut'.:  were  rollected 

for  phosphate,  nitrat  .  silicate,  productivity  m :  thlorophyll  .  m  ire- 
ments.  [h  bservati  for  the  -on  January  to  J  n  r  1\ 

: Barnes  and  Lov  (19<  • ) .  The  information  gathered  n  these  :ruis<  .  was 
;  m  nt<  i  bj  lata  from  'is  Ore)  n  Stal  Ui  Lv  rs  Ity 'i  r  search  vess 
R.  V.  K 1 .  ka  and  A  -ona .  the  Canadian  Nava .  Vessels  C.N.A.V.  ft.  Anth ony  u  1 
Wr.  itt  iiron  t  operating  with  the  Pacific  Naval  Laboratories,  and  the  Pacific 
Oceanograpl  i  •  Group,  and  the  M.  V .  John  N.  Cota-  of  the  U.S.  Fish  and  Wild¬ 
life  Service.  A  summary  of  the  R.  V.  Brown  Bear '  s  cruises  is-  preset:*  ed  it. 
Table  1. 


TABLE  1 


SUMMARY  OF  AND  1962 

CRUISES  OF  R. 

V.  BROWN 

REAR  IN  THE 

COLUMBIA 

RIVER  EFFLUENT 

AREA 

Cruise 

Sal ini t  y 

No.  Period 

Days- 

Mi  les 

Stat ions 

Samp  lot' 

o>-l 


275 

10-27  Jan 

18 

1685 

67 

958 

280 

7-24  Mar 

17 

1^85 

74 

987 

282 

3-  7  Apr 

4 

765 

6 

80 

287 

8-24  May 

16 

1762 

S2 

1093 

288 

9- 19  Jun 

10 

9t>0 

47 

64 r- 

290 

6-25  Jul 

19 

2300 

77 

1134 

291 

2b  Jul  -  13  Aug 

16 

1700 

28 

521 

292 

14-20  Aug 

6 

700 

9 

147 

293 

14  Sep  -  20  Oct 

36 

4155 

131 

2872 

297 

28  Nov  -  18  Dec 

20 

2122 

73 

1613 

1962 

299 

23  Jan  -  7  Feb 

15 

1644 

97 

1500 

304 

27  Mar  -  12  Apr 

16 

1775 

89 

1693 

308 

7- 19  Jun 

13 

959 

79 

1277 

309 

20  Jun  -  9  Jul 

19 

702 

89 

1108 

310 

10-23  Jul 

13 

1506 

31 

647 

311 

24  Jul  -  14  Aug 

21 

2300 

15 

242 

312 

14  Sep  -  9  Oct 

25 

3840 

158 

2887 

’OTALS 

> 

284 

30460 

1152 

19413 

p"  •  ■  •  ■■  Ins  •  '  inn  t 

Convent  iona :  muppi:ig  and  .  ciontific  inv  ■.  t  igations  of  rivers  have  In 
*->>e  pas’  onfii  .1  to  hydrological  studies  of  the  drainage  ;  ystem  oi. 

■  tl  tt  ■  .  •  1  ffluent  :hara  :terlsl  tl 

I  '  ;  •  However,  i  rec  t  years,  nterest  tl  prol  m  ■  n  i 

with  the  lisposal  of  radioactive  wu:  te  and  .  -wage  in  the  op-.-:,  seas  has  give 
rape*  tif  ■  . : :  :  ■  ■  ;  :riptit/  nd  -  retica  L  (  .  . .  HAS- 

® '  •  -  ■  •  :  -  ;  W  lichu  1963).  Tl  lescript  \ 

roit  i  r  . : . :  .  tl  follovin 

is  tics; 

( • )  Ph  ■  fr  h  river  water,  after  entering  tl  sea,  i  tonfined 
trfac  ■  •  "  •  that  ..  met  rs  t  \  .  .. 

nta  listributi  .  isua  j  shaj  f  ■.  p .  um  , 

pr  .  ..  ver  ar  area.  Chau  (  >)  showed  that  Pearl  River 

vat  r  tontributed  as  mush  a.'  2C  percent  of  the  water  in  the  upp  r 
-n  r  yer  Soutl  Ihina  Sea  ]  km  fr<  m  th  Ri\  :  . 

3ate:-  (.933)  found  that  the  Mississippi  River  flowed  from  South- 
t  Pa  Ik  two-d  nsiona  Jet.  rhe  distinguishat  fflu- 

nt  wa  nfined  1  th<  p] . s,  8  kn  fron  tl  ■  th. 

I  .■  v  ■  )  found  tr  f  Mississi pp  wa  t  r  Li  th  surfac 

•  ayr  7^0  to  -V$  km  from  the  delta. 

(•'■)  Tlie  horizontal  distribution  of  the  effluent  water  at  sea 
•1  :  n  1  overned  by  a  mbinati on  of  th  c  ffi  hor  - i r  ;u la- 
tioi  ind  tl  :a  previ  1  wind.  Bol'shakov  (  )  .  .  *ved 

’•he  :  ir-  ct  ion  of  the  plume  of  water  from  the  Dnieper -Bug  estuary 
'  ■ :  n  ■•  ■  th  wind  lir  :tion.  Ph  lireotion  tak  n  bj  th 
.  outhea.  t  nrar.sh  of  the  Kuri  River  was  found  to  change  from  south 
'  rl  ast  ..  Ithi  1  w  lays,  ippa rentlj  in  resj  ns  o  a  wind 

fl  f . .  .  .  . ■  ■  .  :  (D<  br<  vo  'skii,  t  al. ,  L96  ). 

Ssrutor.  and  Moore  (  .033)  found  a  significant  correlation  1.  tweer. 
the  Miss,  issippi  River  plume  srientation  and  wind  c.  irection  win:. 

' : .  rurrenti  u  r  w  and  po radio .  Phej  tone  ided  that  t  he 
wind  stress  was  of  primary  importance  ir.  distributing  the  plume 
lurin  .  peri  1: .  Henrj  ind  Elder  (1958)  found  botl  th  wind 
and  prevailing  -urn.  at  to  have  important  effects  on  the  Mississippi 
plume  orientation. 

(3)  Although  Iehiye  (l-kcO)  concluded  that  river  drainage  has 
r-Lativ  iy  email  influence  on  the  dynamics  and  water  budget  of 
the  adjacent  ocean,  there  are  certain  marked  influences  on  sedi¬ 
mentation  and  biological  properties  directly  attributable  to 
river  ‘ffluent.  Bates  ( i)c>3 )  has  documented  the  role  played  by 
the  Mississippi  River  in  .  -dimentatior.  near  its  mouth.  Riley 
(1937)  reported  a  sizeable  increase  in  phosphate  and  an  associated 
Increai  p  nt  biomass  within  surface  waters  off  the  Mississippi. 

Tliis  report  introduces  the  results  of  physical  oceanographic  in- 
v-s t  igat ions  of  the  Columbia  River  effluent  dispersion  during  ■  and  oY  ' 


hydro!  .  .  ..  ther  background  'eai  wal  r  charaeterlsl 

the 

Quant  it at iv-  -u-timut  at- •  made  of  th-  major  phy:  i  •••»  1  p  i-  r .  'nv  :  •  t  if:--  -t 
the  diaper  ,  and  the  -  dm tnti  :  i  per  n  predict  ten  - 

i  . 


II.  HYDH  1DGY 


PI  wai  1  of  i  Q-  n  .  -n  '  Lum!  RI\ 

pr  .  .  iuare  statute  tulles), 

•  .  percent  r  tt  .  a  (  ,0(  squar  mi  ,  w 

United  State)  1  istitut  it  ut  per  rent  f  1 1  nation1) 
lrainage  ( B  i -  i  .  f  R  una tion,  •••  r  ■  . 

f  Was)  .  nd  western  Mot 

Nevad  l  Utah  (Fig.  1 ) .  Tl  ir  the  ri\ 

f  British  C  imbia ,  meters  (2, <  f  t  abo\ 

w.  ••  ■  :cr:  (•••  1  :r: i  .  1  i::v.v'h  d'C:?:.!  .  •  d,  -if.  ;•  :  ' 

•  ••  :  -orr  ••  th  fat--  W - 1 .  1 1  i : if :  ,  •  •  •  .  r.virr  ' 

the  Wasl  n-Or  mdary  wher 

'  P  ••  at  . 


*r  com- 
I  Lght  y- 
i  the 
.  The 
Idah  . 

W  n  ■ 

•  is  i ■  a 
. 

Lag  th- 
h  th 


rhe  River 1  ar  t  tr  ■  R . id  . 

.  r  r  ■■  ’  •  Ri\ . 

H  mette  .....  .  t ter  tv  rlbut.ar  •»  v 

f  The  Sail  .  1 1  point  al  wh id  I  C  mbia  penetr  th 

V.  ■  •  I  .  :  :"U\-  . 

he  am  mean  die charg  f  1 1  '  ml  River 

,  ' .  mp  stations  re i  *t  r  nd 

sec”  timated  by  Leoj  l  (1962).  rh  is  rut  ‘f  oi  r 

::  :  acre-f  t  j  r  j  . epres  nt)  ’7  I . ent  >f  the  tot  .  .  3  nd  in- 

gaug  d  freshwater  dir  -harged  into  the  Pacific  Ocean  south  of  th  Strait  of 
Fuc  ind  includii  l  R  River  at  • .  JO'N  at  *  bout  6( 
percent  of  th  tota  ;auged  md  tug  1  lii  :hari  tc  bu  including 

ff .  ,  :  ■  f  S  : . Lsc<  Bay.  rhe)  percent)  ir  wt  r  :on- 

.  ■  yearly  period.  Difference)  •  soi 

be  tv  peak  rui  ta  ri\  ind  that  f  lo  i  River 

.  •  1 .  rtuations  th  r  ntributions  of  hwat  r. 

f  tl  '  unbia  R  . un  M  '  Jun  re- 

■  ■  now  n  ■  .  wl  ■  1  for  tl 

the  R  R  .  ■  ir  luring  th  /y  pr  :i]  tatioi  j  ri od  r  n 

;  . .  February.  I  ir  period  s  one  f  n  mun 

nd  part  f  tl  w  r  isin  f  tl  Co  uni  R  .  r.  Durii 

winter  month;  ,  the  Columbia  River  rontributes  about  60  percent  of  the  fresh 
water  Ided  t  -  I  :ean  ■  ;  Lon,  but  dur  tl  .  uam  r- 

ear Ly  fa  per  l,  wh  in  th  riven  r  .  ;ontribute) 

mu  cl  .  .  Phe  . nofl  fron  son  'ivers 

ft  r  ■  .  ..  ...... ..  ..  ..  •  •  -  un,  Durii  t  ■  wint  r 

River  lei  nix  w  [tl  idmal  v  rivera 

.....  .  ....  liffi  cult  to  distil  i  r  I  ntity. 


BRITISH  COLUMBIA 


CAPE  FLATTERY 


WASHINGTON 


•  cob  A 

Sf  AMOUNT, 


IDAHO 


OREGON 


CAPE~MEfcDOCIN(r 


2500 


NEVADA 


PT  CONCEPTION 

I  W. 


in  Mini  udjaoont  sen. 


cT'  'SOO  fm 

2000  fm.  >■ - 


PRAINAGI  TO 

SAN  FRANCISCO  BAY 


CALIFORNIA 


J  how  "  'omputed  t  i  r  for  the  n 

rairuiy  •  iw  ,-n  t  hot-  riv.-v  fr.  •:«  th.  fd  r.ait  v  J  .  • ,  :  F.  .  to 

"in  ’  ■  '■  Bii.v.  Th.  V.M  hi  :  h> -w:  u  re  •  ,  ,t  i  ;‘r  : 


lotal  hi.  ehaty.e  Gaiii •  .1  1) 1  'bar 


(imHua.,) 

O'  HI’A'  T  j 


v  re  Pu  ami  -  reoipit 

reSP  ;  '  ;  |,;  •  »  :  DA  gauged  d 

reaped  i\  .  IYi  wer  n  from  ; 

U.S,  Army  Cor  Stigineers. 
about  {'•  per  a' tit  of  th"  tola  ur 
from  the  10-  to  to- war  -,v. 


G1,  II,  d  d  i  >ha  r,  "  f  i  I, . t;s  . 


"  ‘  1  '!  01-".  " 

' , 

A  tv  .  1 1  • ;  L*.  .  »•,  .  • 

•  Ian  i m  -ter 

,  .  .  :  - 

1  :Sur\  ,v  map:  . 

TABLE  S 

* '  ■  *  » *  "  t 

MEAN  ANNUAL 

DISCHARGE  OF  COASTAL. 

RIVERS 

Area  (ini'  ) 

Dir. char  o 

Ru  i-.ef  f 

Strait  of  Juan  .1  •  Put*  a 

(ft 

( :'n  -ti  ■■  ) 

t  o 

Gray  Harbor 

' .  a  \ 

20,t'ii' 

11 ! 

Gray.  Hart-  -r 

.  1  -  V  - 

l .  ’ ,  V  1 

t-0 

Wi  a  pa  Ba  > 

oi(S 

h.Oi'a 

73 

00U1MB1A  RIVER 

S  o.tX’O 

Y-i'.OOO 

39 

Fr.  :-i  C.  unu  i  K  i  \  r 

(b- 

1  ow  Th.-  D- 

to 

Umpqua  lii\  r  Bar  in 

.  ,  : 

.  '0 ,  ■*  . 

t, 

t  “4 

Ur pqua  Rlv-r  Hr  in 

Fr.-  Umpqir  Ba  i  u 

■t .  s- rf 

30 

R<  ■  -  i  I  n  B 

\  iid 

i-0 

Rogu.  •  and  rilinoi.  R  \ 

Fr  -m  Ro.o.ie-Il  .  '  no  i  r  Bar  in. 

'  .  i  >  i 

e  ,  o 

3>' 

to 

Y .  YT‘ ' 

vr.'t.s 

so 

,  ,  Til"  data  .in -si  to  compute  the  percentage  •.  .at  ri  Initio-  i\ 
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Geohv.i.-il  Wat  -r  Dap-r.  .  IV- low  Th.  Dalle,  the  Columbia  i  not  n,,f-  i  and  only 


■  i 

i  i 


k 
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K  tkital  ri\  rs.  The  torreotion  factori  a  no.  qua!  d  ■  •  tomput 

daily  total  d  isoharp'  of  the  Columbia  R  i  ver  at  A. -tor .  a  .ar>  pro.  at  .• 
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Dlre.-t  me.a  suromen  t :  .at  h  ■  mouth  of  tlv  Columbia  Itivor  for  part:  of  Hay, 
Juno,  and  -pt  ember  (a  total  of  a'  ,ia,\.  )  are  a\.a  i  Cab  i  •  for  >>•,•>  (ll.f .  Aam\ 
Eng  i  neors ,  1'H'O).  A  owtwlnoit  of  th.-  daily  d  i  soh.aiy.o  oompiutod  a.  do.orib, 
abov  •  with  that  measured  t-y  the  U.C.  Army  t’orp.  of  Kng  i  n  vrs  how.  ood 
agr  'enietit  .  If  a  ‘.-day  la,,  period  i.  applied  t  ,•  m  •.asuromotit  .  :•  ad  •  at  th 
•  t  available  .gauging  st.at  lotu  .  th.  ootn pitted  daily  dlsohaty.  •  it  tb  mouth 

within  i  percent  of  that  r  |  iri  d  by  \rmy  Ei  _  . i  .  Tit 

’•  it  •  tomp&r i s on _ vnr i ed  between  •  .c  •  and  16.7  x  : c '  t:.  '  .  to  t,  . 

and  ‘’00,000  ft  see  .  respent  Ively),  The  Dalles  Is  w,'  km  frotp  1 1 1  river 
moutti  and  the  -day  lag  would  suggest  a  it.  an  flew  of  1 1  em  .  «e  '  ;  how  \v. 
the  .apparent  b-.iny  lag  !:'•  in  part  an  arttf.aot  ore, a  ted  by  the  loelmiqu  of 
applying  •orr>,ot  ion  faotors  for  the  lower  basin.  Re  l  inu.v  may  b  plan  d 
the  daily  dlsnhfirge  It’  a  '-day  i.ag  is  .applied;  hewovop,  the  diff.ovr.o 
Ins  igtii  f  I  oanl  for  our  purpose  .  and  till:  lag  Is  not  i  noorporat  od  in  the  da 
dtsnh.urge  tabu  In  t- ions  shown  in  Figure  '. 


>  I 


Tli'-  results  of  the  d  i  i  charge  oomputat  ion:  l’or  i  'cl  and  J<  nr  h  w:i 
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Rogu  •  Ri  v  a*  in  i"  vi .  000  km'  ffV'i.OOO  mi'  ).  Tho  compu  t  ■  *«i  annua  1  ih.'h.ur 
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Tho  rolo  of  ground  water  in  holiovod  to  bo  of  Little  importance  in 
the  quantity  of  laily  lit  tharge j  however,  the  Level  fluctuation,  ii  ol\  d- 
.  o  1  idr  compos  i  t  ion ,  and  circulation  of  ground  valor  an*  s  Ign.i  f  i  ••ant  v.uri- 
ibl":  a  f  foot  ittg  tho  quality  of  the  w.utor  along  th"  reachoe  of  t  In  R :  v  ■ " . 

part  ieularly  in  th  •  lowland  aiva.  .  From  tho  oceanographic  rtan.ip.  in*.  ach¬ 
ing  by  ground  wator  may  contribute  trace  subst ances  no. •or;  ary  for  th  intr- 

mi 'diary  mo  t  a  ho  1  i  run  of  plant;'  and  animal:-.  Tho  rolativoly  low  productivity 
-luring  tho  wintor  rapid  runoff  :  •urou  may  bo  duo .  in  part,  to  a  ..-w  -r 
oonoont  ration  in  wintor  than  in  rummer  of  truce  rubrtnncor  which  cat  a  iy. 
growth  and  production.  Tho  corn*  la  lion  hotwoon  fror.  hw.at  or  and  high  produc¬ 
tivity  in  tho  con  ir  d  ircurrod  by  Attdorr-on  (  1  *>('•': )  in  a  compan  ion  \ .  m;  of 
thl:  report. 


The  Washington  (State)  Pollution  Control  Commirr  ion .  th*'  Poparttmu  : 
of  Conrorvut  ion  and  tho  U.G.  Geological  Survey  (  1  * H •  1  )  Joint  ly  report  tho 
moan  discharge  of  dissolved  solids  at  The  Da  l  ion  to  bo  0l»  ppm.  Th  -  two 
major  components  of  the  dirrolvod  material  are  hi  carbonate;'. .  comprising 
about  '*  porcont  of  tho  total,  and  calcium,  making  up  about  !•'»..<  porconi . 

IV  enure  the  ratio  of  dissolved  conr- 1  .i  tuent  r  in  the  river  differ  from  there 
in  tho  con.  tho  apparent  value  of  salinity  determined  from  conductivity 
will  be  inaccurate,  but  this  error  is  negligible  except  in  urear  of  low 
salinity  along  tho  lower  river  course.  Tho  month- to-month  changes  in  con¬ 
cent  rat  ion  of  the  dissolved  solids  are  out  of  phase  with  the  river  stage; 
i.  ..  ttie  higher  ooneontrnt  ions  of  total  dissolved  sol  idr-  occur  when  the 
discharge  is  low,  but  the  ratios,  of  the  individual  dissolved  solids  appear 
to  remain  rather  constant  throughout  the  year.  Nit  rate  ;  .  vice  as  abundant 

in  the  river  water  during  the  into  winter  and  early  spring  than  it  : :  ir. 
summer  and  fall.  A  discussion  of  the  relation  between  river  water  nutrients, 
and  ambient  ocean  water  nutrients,  is  given  by  Stefans  son  and  Richard:  ( ’ek-j' 
and  the  analys.es.  of  nutrients  and  other  cons  t  1 1  went  s  in  the  River  .are  given 
by  :'.y  Ives  ter  and  Carlson  (idol). 


III.  WEATHER 


Precipj  tat,  ion  and  Evaporation 

Thi'  mean  annual  precipitation  in  the  coastal  and  immediate  offshore 
areas  of  the  northeast  Pacific  Ocean  influenced  by  the  Columbia  River 
effluent  Is  istimated  to  exceed  L20  an  yr  Most  of  th  precipitation 
falls  during  the  winter  season,  the  summer  months  being  generally  dry. 
Precipitation  exceeds  evaporation  and  essentially  the  entire  ocean  area 
under  consideration  is  me  of  net  dilution.  The  surplus  of  precipitation 
is  least  in  the  west  and  south,  and  increases  to  a  maximum  of  more  than  CJ0 
cm  yr  '  off  the  coast  of  British  Columbia  and  southeastern  Alaska.  South 
of  the  Columbia  River  th--  surplus  decreases  to  near  zero  at  San  Francisco 
Bay  (Jacobs,  19!?l).  The  precipitation  over  the  Columbia  River  irainnge 
ba.  in  and  contiguous  coastal  drainage  bat  Lm  i:  a.-  variable  a:  the  topography. 
In  the  Columbia  River  Basin  alone,  the  yearly  average  varies  from  less  than 
20  ■  over  300  an,  depending  on  Location.  The  mountainous  areas  receive 

generally  heavy  precipitation,  much  of  it  in  the  form  of  snow,  whereas  the 
plateaus  lying  in  the  rain  shadow  of  the  Cascade  Mountains  receive  little 
rain  or  snow. 

On  the  average,  during  the  summer  the  surplus  precipitation  at  sea 
is  about  0.1  cm  day  and  during  the  winter,  about  0.3  cm  day  ,  according 
to  estimates  by  Jacobs  (1951). 

The  surface  exchange  of  water  by  precipitation  and  evaporation  in 
the  Columbia  River  effluent  area  should  differ  but  little  from  that  of  the 
.-unbient  ocean  water  and  -an  probably  be  neglected  in  the  first  approximation 
of  the  budget  of  river  water  at  sea. 


W  i  nd 


Winds  on  a  regional  basis  affect  the  general  water  circulation  in  an 
area  and,  on  a  local  basis,  both  the  movement  and  mixing  of  surface  waters. 
The  seasonal  cycle  of  winds  over  the  northeast  Pacific  Ocean  is  largely  de¬ 
termined  by  the  circulation  about  the  North  Pacific  high  pressure  area  and 
the  Aleutian  low  pressure  area.  During  the  summer  months  the  high  reaches 
its  greatest  development.  In  July,  the  center  of  highest  pressure  (approxi¬ 
mately  1,02!?  millibars)  is  located  near  30°N,  150°W;  and  an  average  pressure 
in  excess  of  1,01!?  millibars  prevails  over  most  of  the  northeast  Pacific- 
Ocean.  During  this  period,  the  Aleutian  low  is  almost  nonexistent.  This 
pressure  distribution  favors  predominantly  northwest,  and  north  winds  over 
the  coastal  and  near  offshore  areas  of  Oregon  and  Washington. 

In  October,  the  Pacific  high  extends  from  the  United  States  coast 
across  the  Pacific  Ocean  and  onto  the  Asiatic  Continent,  reaching  a  maximum 
of  1,020  millibars  pressure  In  the  vicinity  of  j0°-3!?oN  and  13!?°  -  lU0°W.  The 
high  weakens  with  the  approach  of  the  winter  season  and  by  November  it  i. 
little  more  than  a  weak  belt  of  high  pressure  lying  between  the  Aleutian  low 
and  the  equatorial  belt  of  low  pressure.  The  Aleutian  low,  with  pressures 
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lower  than  1,002  millibars,  appears  on  the  charts  as  a  permanent  system 
located  near  50° N,  l80°W,  but  in  reality  is  a  system  of  migratory  lows  fre¬ 
quently  moving  through  the  area  during  the  winter  season.  These  traveling 
depressions  moving  eastward  cause  considerable  day-to-day  variation  in 
pressure,  particularly  in  the  area  north  of  40°N.  The  winter  winds  are 
frequently  of  gale  force  (5  to  8  percent  of  the  observations)  and  range  in 
direction  from  southeast  at  the  coast  to  southwest  in  the  offshore  region. 


IV.  ESTUARY  AND  NEARSHORE  ENVIRONMENT 


The  initial  mixing  of  the  river  water  with  seawater  occurs  in  the 
estuary  and  nearshore  environment  where  tidal,  surf,  and  hydraulic  processes 
prevail.  The  relative  importance  of  the  mixing  sedimentation,  and  bio¬ 
logical  processes  in  the  estuary  and  nearshore  environments  were  not  studied 
in  1961  and  1962  by  the  University  of  Washington.  Model  studies  are  being 
conducted  by  the  U.S.  Army  Corps  of  Engineers  and  the  U.S.  Geological  Sur¬ 
vey  at  Vicksburg,  Mississippi.  These  studies  include  movements  of  dissolved 
solids  and  bed  load  in  the  estuary. 

A  study  of  the  salinity,  temperature,  and  velocity  of  the  Columbia 
River  between  Beaver  Army  Terminal  (mile  52)  and  the  River  mouth  was  made 
by  the  North  Pacific  Division  of  the  U.S.  Army  Corps  of  Engineers  (i960). 

Data  were  collected  in  1959  from  a  number  of  cross  sections  at  different 
depths  and  under  conditions  of  low,  medium  and  high  river  flow.  According 
to  these  measurements  the  River  hydraulic  current  is  reversed  on  flood  tide 
during  all  seasons  of  the  year  up  to  40  km  from  the  mouth.  During  the 
September  low  discharge  period  of  1959>  flood  tide  reversed  the  River  flow 
at  the  highest  upriver  measuring  section,  8k  km  (mile  52).  The  intrusion 
of  salt  water  extended  upriver  past  21  km  (Tongue  Point)  but  less  than  32  km 
during  September  1959  and  not  past  19  km  (Astoria,  Oregon)  during  the  high 
discharge  period  of  May  and  June. 

The  Pacific  Ocean  tides  can  be  measured  upstream  as  far  as  Bonneville 
Dam  where  the  diurnal  range  is  0.2  m.  The  tidal  prism  volume  for  the 
Columbia  River  from  the  mouth  to  Bonneville  Dam  (23^  km  upriver)  is  about 
O.85  km  .  Columbia  River  water  enters  the  nearshore  environment  in  12-hour 
pulses  as  clouds  of  low  salinity  effluent. 

In  1963,  examination  was  begun  of  physical  processes  and  the  scale 
of  phenomena  in  the  nearshore  environment.  The  observations  indicate  large 
tidal  influence  and  an  initial  mixing  in  the  estuary  of  about  one  part  sanw 
water  with  two  parts  water.  The  currents  near  the  mouth  of  the  river 

directly  affect  the  local  transport  and  mixing.  At  the  mouth,  midway  be¬ 
tween  north  and  s<j>uth  jetties,  the  measured  currents  near  the  surface  were 
about  300  cm  sec  seaward  on  ebb  and  120  cm  sec  upriver  on  flood  during 
June  and  240  cm  sec”  seaward  on  ebb  and  180  cm  sec  upriver  on  flood  in 
September  (U.S.  Army  Corps  of  Engineers,  i960).  From  5  years  of  discontinu¬ 
ous  measurements  (1915  to  1920,  with  1915?  1919?  and  1920  nearly  complete) 
of  surface  current  by  drift  pole  (Marmer,  1926)  it  was  possible 

to  resolve  the  hourly  observations  into  tidal  and  nontidal  components  at  the 


. 


1? 


Columbia  River  Lightship  located  tan  southwest  of  the  south  Jetty.  Tin- 
tidal  curr  nt  is  -lookviee  rotary,  semidiurnal  witii  little  diurnal  inequality. 
Strength  f  flood  i .  •  cn  sec  )  come;  about  one  hour  befor<  Low  water.  Flood 
sets  north  aster ly  and  :•  s<  ithwesterly.  The  nontldal  turrent  setting 
westerly  and  southerly  ranged  from  a  monthly  average  of  15  cm  cec  in  March 
to  3b  cm  dec  ‘  in  June.  There  is  some  correlation  between  discharge  und 
speed  of  current  at  the  lightship.  The  yearly  average  ijet  was  west  and  the 
average  of  all  measurements  (over  20,000)  is  19  cm  sec  .  Hydraulic  cur¬ 
rent  s  are  considered  by  Taka  no  (195*0  to  diminish  by  an  order  cl-  magnitude 
a:  a  distance  seaward  equal  to  several  tines  the  width  of  a  riv  r  mouth. 

This  concept  is  supported  by  these  measurements  off  the  Columbia,  and  hence 
beyond  20-30  km  from  the  river  mouth  geostrophic  currents,  wind  stress, 
diffusion  and  biological  transport  art-  probably  the  most  important  processes 
aff -sting  dispersion  of  river  water  and  river  borne  materials. 


V.  PHYSICAL  CHARACTERISTICS  OF  THE  OCEAN  REGIME 


Ambient  Water  Structure 


The  northeast  Pacific  Ocean  is  a  region  of  net  dilution,  in  which 
the  precipitation  on  the  sea  surface  and  adjacent  land  masses  exceeds  the 
evaporation.  The  land  runoff  acts  as  a  line  source  of  freshwater  extending 
from  southern  Alaska  to  the  coast  of  northern  California.  Along  this  line 
the  Columbia  River  is  by  far  the  largest  single  contributor. 

Tin-  area  of  the  '.Northeast  Pacific  Ocean  subject  to  t lie  greatest 
influence  by  the  Columbia  River  outfall  is  bounded  by  latitudes  -t0°  and 
5 cf  H  from  the  v..  t  ( L24°W)  west  to  :  \’°w,  and  is  subarctic  in  character. 

Rev  lews,  of  the  oceanography  of  the  subarctic  region  have  been  given  by  var¬ 
ious  authors  and  most  recently  by  Dodimead,  Favorite,  and  Hirano  (l'X'3)  and 
Uda  (l9t'3)»  It'  the  upper  kilometer,  the  vertical  thermohaline  structure  is 
•hornet er iced  by  three  layers  or  zones  (cf.  Doe  1955.  Dodimead  195''. 

Fleming  1031').  The  first  of  these  zones,  the  surface  zone  of  75-100  meters 
thickness,  is-  somewhat  diluted  because  of  the  runoff  and  net  precipitation; 
the  dilution  is  not,  however,  accompanied  by  large  changes  of  salinity  with 
depth  in  the  .'.one.  The  temperature  within  the  layer  decreases  markedly 
with  depth  in  summer,  but  in  winter  is  nearly  uniform.  The  second  zone, 
a  quasi-isothermal  halocline,  shows  a  marked  increase  of  salinity  with  the 
depth,  and  falls  approximately  within  the  100-200  meter  depth  increment. 

In  the  deeper  waters,  which  comprise  the  third  zone,  the  temperature  slowly 
decreases  and  the  salinity  slowly  increases  with  depth. 

Tully  and  other  investigators  (Dot'  1953,  Tally  and  Dodimead  1957, 
Dodimead  19s''.  Tully  and  Barber  10t'0)  have  shown  that  the  salinity  at  the 
bottom  of  the  halocline  is  very  nearly  constant  (330  +,  0.1  '  /  oo )  throughout 
the  year.  The  values  at  the  upper  boundary  of  the  halocline  (which  varies 
i tv  depth  from  75  m  to  150  m)  change  throughout  the  year  0.2  /'oo  or  less  at 
any  given  loeation.  In  tin'  area  affected  by  Columbia  River  discharge  the 
salinity  at  the  upper  halocline  boundary  is  approximately  32. 7  /oo,  but  in- 
■reas.es.  toward  the  west  and  decreases-  toward  the  north  (Figure  3). 


(  0  mile:  per  lay ) ,  .  ettitig  south  i  r.  .  pring  fiiul  .  unit.-  Dr:!'* 

:  tud  ies  L_\  Thompson  anti  Van  Cleve  ),  Dod.n-id  and  H  .  .  !  st  r  ( 

Ren:  and  Wyatt  and  Kujala  (’•>(.')  •  evoior:  •  ...  .  •  ..  ,  . 

■  1  •  ■  iy nam ic  studies .  Ih<  ir 1ft  bottle  stud let  f  TL  mj  -  1  Van  C lev 

v  >36)  and  iynamic  :a  tulationi  L  \  fully  ( 1938 >  si  ■«  1  w  which 

was  northerly  in  the  .  pring  and  south  -rly  in  the  rummer,  probably  :•••  lat  d 
to  th  lasonal  :hange  In  the  wind  pattern.  Dodimead  (  -  )  hai  liscu  I 
th>  year-to-year  variations  in.  the  geos  trophic  -irculation  from  l'» -\>  t. 

: 


Ir.  t in-  prerent  study  a  dynamical  analysis  tuu'  beer,  made  of  vtie  wnt-r 
•ni.bi  •  •••.;  to  th  ■  ’oa:  :  r  ef  Wa:di  Ittf-tt  on .  Ore  (.'on,  and  Northern.  California  and  th 

.  rapt  art  shown  it  Figures  ‘  to  l4,  Geostrophic  tu . 

(•tie  vni^er  movement:-  associated  with  tin1  internal  distrinut  ion  of  mass) 
wer  .  luced  fron  :u lat ions  of  topography  of  the  sea.  Th  L 000-decibar 

.  irf  vai  :  u  the  ••  ference  Lev*  1  (level  of  as*  um< *d  no  mol  Ion) . 

Ca i ru iations-  were  made  for  all  Browr.  Pour  rrui.es  during  19*1  and  1  >* -cx- 
•  pt  Crui:  (geology  -mire)  and  Cruise  309  (biology  •ruls-el.  The 

t  opographies  shown  are  deviation;,  from  the  geo  Id  surface  and  ne<  1  -,-ct  the 
ff  ■  ■  of  atmospheric  pre:  ire .  These  charts  repres  nt  tl  first  season- to 
lyt  m  ;<  nalysii  f  the  offshon  area:  adjacent  1 .  the  Mash ii  ton, 
and  ,  ■  C  Lifornia  toasts.  Throughout  this  analysis  standard 

nn  .  :  .  .  f  mtputit  Ij i  am Lc  h<  ht  an*  ma  1  >s  were  us  td .  Va  u  is  for  1 1 

shallow  in.  hor  station,  were  -omputed  by  a  method  of  extrapolation  (He .  lat. 
Hansen  193M  or  by  the  technique  of  plotting  the  specific  volume  anomaly 
ri  u  the  Logarlthn  f  lepth,  u  1  .  :rit  d  by  Bennett  (1959)*  or  by  both 
methods.  H  stations  leeper  than  ,000  meter:  are  circled  on  the  charts 
Station:  f  .  .  than  JOC  1  iter:  lepth  w* ire  not  used  in  the  extrapolation. 


Tiie  cont  our  intervn.,  .'  dynamic  oent  imet  ors,  used  in  ,vnts:;rs.s  tin 
topographic  charts  may  be  open  10  criticism  ir.  view  of  the  probable  'from 
1 ;  measurem  nt:  at  th*  '  c,  .••  imploy  -a  (cf.  Reid  .0  •>:.  Lpp  r  95  >)« 

Large  f luctuat  lor.s-  can  also  occur  in  individual  dynamic  height  calculations 
because  of  local  short-term  fluctuations  in  the  distribution  of  mass.  ir.  the 
wat  r  loluan.  it  ha:  beet  found  that  these  f  luc: uat  i ons  are  oft--:  of  the 
same  magnit  td  as  tl  :ot  tour  interval.  The  probable  error  of  *  2  dyn.  cm 
is  based  on  repetitious  sampling  during  these  observations  and  The  results 
repetitious-  measurement  s-  at  station  PAPA  (Fofonoff,  personal  communion!  ion) 
The  short-term  fluctuations  were  probably  due  to  the  passage  of  intense 
storms  and  internal  why  -.  as  well  as  errors  in  measurement.  Defant  (lose), 
it:  examining  a  series  of  geopotential  topographic  chart  s  prepared  by  Seripp 
Institute  of  Oceanography  for  an  area  off  the  coast  of  California,  found 
that  many  of  the  highs  and  lows  could  be  attributed  to  internal  waves. 
Wooster  and  Taft  (ld‘>8)  have  shown  that  a  probable  error  of  *  1  dyn.  cm  can 
be  expected  considering  the  precision  of  measurement  only. 


Tiie  reference  level  of  1 , 000-decibars  was  chosen  because  of  tin' 
limited  number  of  observations  in  deeper  water.  Reid  (i9ol),  Dodimead  (  oV 
and  Bennett  (1959)  have  suggested  that  1 ,000-decibar  level  is-  adequate  for 
show  the  lirection  of  flow  in  the  eastern  subarctic.  It  is  Improbable 
that  a  uniform  reference  level  exist:  in  a  north-south  direction  in  the  sub 
arctic.  Neumann  (l 955)  has  shown  that  the  reference  level  rise:  from  the 


Equal  ••  •  ov&r  l  th<  Pole!  In  the  Atlantic  Ocean.  Dr  ■  measurement  mad 
by  the  Brow:.  Bar  in  19trl  in  the  effl.:- ■  nt  area  indicate  a  mrr  :.t  at  i,CCC 
n  '  rs  f  2  .  In  view  f  thii  lit  nil  si  n,  :harti  f  Lyru  m 

topography  0  over  1,000  decibars  ar-  probably  slightly  in  -rror  in  som 
parti  f  tl  area;  but  Li  felt  that  they  adequately  lepict  tl  rosi 
features  f  t hi  speed  and  direction  of  water  movement.  ..  onstant 
decibar  reference  l>vl  is  inadequate  for  computations  of  mas:  transport  in 
this  area. 

During  January  1 ,  a  current  wa:  observed  to  s<-t  northward  at  ap- 
pr  ximat  Ly  sa  sec  '  (Fig.  ■  ) .  This  flow  intern  if  led  t<  gr  at  than 

:m  s  -  ’  during  March  (Fig.  •  ),  and  reversed  iirection  luring  May  (Fig.  T). 
The  radical  :hange  in  :urrent  lei  Li  associated  with  a  shift  the  predomi¬ 
nate  wind  from  southerly  to ^northerly  and  during  the  summer  a  southerly 
turret  t  av  (raging  as  sec  '  pr  1  >minat  «i.  The  cure-  nt  pat  t  rn  :hange  fr  m 

winter  to  summer  ii  :omplicat  d  by  jyclonic  and  anticyclonic  Idles.  Then 
is  a  suggestion  of  a  cyclonic  eddy  in  the  topography  of  Figure  6  and  this 
Idy  appears  bett  r  developed  in  Figure  J ,  The  anticyclonic  Idy  or  high 
near  -•  shown  in  Figures  7  to  10  is  well  documented  by  numerous  observa¬ 
tions  and  appears  rather  consistently  throughout  the  summer  and  autumn, 
hence  is  probably  not  an  artifact  associated  with  the  lack  of  synopticity 
Ln  sampling  and  r  internal  waves  (Defant,  L950)*  Reid*  Schwartz  lose,  and 
Brown  ( 1CX  s )  confirmed  t.he  existence  of  a  permanent  eddy  off  the  California 
:oast  ..  Lng  parachute  irogues.  The  Irogue  movements  agree  with  what  might 
be  expected  from  dynamic  topography,  which  suggests  that  the  eddies  ar 

ostrophi  ally  Ln  balance.  Reid  and  co-worker .  ■■  st<  1  that  the  iddies 

pr  bal  y  result  from  the  horizontal  shear  between  two  rurrents.  McEwen 
(lo4b)  studied  the  dynamics  of  similar  eddies  and  concluded  that  the  de- 
rrease  In  v  Locity  and  increase  in  si:  not  d  on  successive  :ruises  are 
probably  successive  stages  of  decay  from  ar.  initial  state  of  maximum  inten¬ 
sity  and  small  extent. 

The  dynamic  topography  derived  from  the  early  fall  cruise  (September- 
Ootober  idb : )  is  shown  in  Figure  10.  Velocities  as  high  as  20  :m  sec  ’  ar 
evident  in  the  cyclonic  eddies. 

The  response  of  the  distribution  of  mass  to  the  autumnal  shift  in 
winds  from  predominantly  north  to  predominantly  south  is  shown  in  the  geo¬ 
potential  topography  for  December  :*v;  (Fir:.  Ll).  Hie  northward  setting 
Davidson  Current  appears  to  be  in  an  early  stage  of  development  but  .he 
isopleths  inshore  are  based  on  Helland-Hansen  ( 193*0  extrapolation  over 
large  distances,  hence  (he  depiction  may  not  be  accurate.  The  circulation 
pattern  of  the  surface  water  is  less  complex  in  mid-winter  as  shown  by 
FJ  jure  12.  Here  a  trough  is  evident  about  L30  to  L85  kilometers  offshore, 
with  a  northerly  rurrent  inshore  and  a  southerly  current  seaward.  The  speeds 
of  the  two  opposing  flows  are  about  equal  at  10  cm  sec  . 

The  southerly  winds  persisted  until  after  the  next  sampling  in  April, 
19o2.  Figure  13  shows  the  northerly  (Davidson)  current  fully  developed 
along  the  -oast.  During  May  and  June.  1062.  the  winds  shifted  from  southerly 
to  northerly  and  the  spring  conditions  again  appeared  with  a  nearshore 
southerly  current  and  an  offshore  northerly  current  (Fig.  Ik).  Thus-,  in 
general,  the  average  surface  current  pattern  along  the  coasts  of  Washington 
and  Oregon  is  northerly  at  10  to  20  cm  sec  ‘  in  the  winter  and  southerly  at 
5  to  20  cm  sec”  in  the  summer. 
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•  BROWN  BEAR  CRUISE  280(7-24  MARCH  1961) 
O  STATIONS  DEEPER  THAN  IOOO  METERS 


Fig.  6- -Geopotential  topography,  °/lOOO 
7-24  March  1961. 
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H 1 g  8—Geopot.ential  topography,  °/ 1000  lieotbara,  Hnvn  Hoar  Cml 
9-11}  June  1901. 
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Fig.  11 --Geopotential  topography,  ®/lOOO  decibars 
28  November-l8  December  1961. 
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STATIONS  Of  f  t'l  H  THAN  1000  Ml  '  I  RS 


Fig.  12 — Geopotential  topography.  °/lOOO  decibars,  Brown  Bear  Cruise  2QQ 
23  January-7  February  1962. 
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VI.  CHARACPERIt'-T ICS  OF  THF?  EFFLUENT  IN  HIE  OCEAN 


Areal  Dint  ribut  ion  of  t  ho  Effluent 


Till  areal  dir-tributionr-  of  tho  River  effluent  for  urUv.t;  t  iv.,. ■;  dur- 
ing  l‘H-1  and  oH\'  are  :--hown  in  Figures  ]*>  through  A.  d  treurr-.\l  above , 

the  l.'.  -  "  oo  ir-ohnline  war-  found  to  bo  a  reasonable  limit  for  tho  deiinea- 
t  ion  of  t.ho  plumo.  Titer  o  ohart  r  of  surface  salinity  show  ,\u  i  t  o  /'.oar'.; 
that  during  tho  rummer  period  of  predominate  uorthoriy  winds  tho  rivor  phi., 
rt  retime-  toward  tho  r-outhwort  ;  and  during:  tho  winter  rear-on  of  prevailing 
r-outhor ly  winds  tiio  piumo  appoarr-  ar-  a  band  of  low  ralinity  water  ad, lucent 
to  tho  Washington  ooart  .  blending  into  tho  freshwater  off  ha"  1  fro::;  tlv 
ooartal  rivorr  and  Strait  of  Juan  do  Fttea.  Definite  limit;  of  tho  p'.amo  ar- 
dororibod  by  a  fixed  salinity  aro  subject  to  r-omo  error  booauro  tho  r-r.rfaro 
water;  -nay  .  coni  ;un  i  nat  >  1  by  ora  pivoipi  t  at  ion  or  '.oca  land  brnlango. 

Tit ir-  ir-  especially  t  rue  during  the  winter  rnonthr.  vlton  the  ooarta  1  rivers 
are  high  and  t  lie  Columbia  ir-  low.  Only  part  of  t  Ho  nearshore  bolt  of  low 
ralinity  water  duritv  tho  winter  ir  duo  to  effluent  from  tho  Columbia. 

Kao torr  r-ueh  ar-  there  r-hould  bo  considered  when  examining  tho  gro:  .  moath- 
to-month  changer  itt  the  plume  ar-  dororibod  below.  Furthermore,  tho  efflu¬ 
ent  cannot  be  detected  after  throe  ordere  of  m^ugnitudo  of  die.;  tior.  b\  roa- 
wntor  booauro  the  limiting  precision  of  tho  r-a  1  inomotorr-  urod  r  dot  ermir  i  ng 


water  booauro  the  limiting  preoir-ion  of 
‘•he  ralinity  ir  about  . COh  '  oo. 


During  the  winter  •miser-  of  l‘Vl  and  1 (Figs.  la,  It-,  .'.'.and  . 

the  Hivor'r  effluent,  oombinod  with  the  local  runoff  from  the  ooartal  rivorr 
extendr-  from  approximately  V  xm  south  of  the  mouth  of  tho  Columbia  to  north 
of  the  Strait  of  Juan  do  Fuoa,  ;r.  a  narrow’  coastal  bolt  about  ‘>0  km  wide. 
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In  the  .  prill#;,  Ma>  1 96 1  and  April  1962  (Figures  17  and  .  ••  ),  the  dis- 
•harge  t'roir.  the  Columbia  is  no  longer  so  -lonely  confined  to  the  coast  by 
prevailing  winter  winds  and  isi  >cia1  d  currents.  During  arlj  .  urd  r,  vfa  i 
■ h  •  astal  runoff  :.  a  minimum  and  the  Columbia  a  maximum,  th  w  salinity 
wat-T  is  readily  identified  with  a  single  source  and  the  affluent  plume  is 
rather  sharp 1>  1  Lneated  by  the  32.5  /  isohaline.  Figure)  and  25 
■learly  show  the  plume  extending  some  390  km  southwest  with  an  . •  I'fshor- • 
boundary  about  210  km  from  the  -oast. 

Just  north  of  the  River  mouth  in  June  and  July  th  ■:*,  is  a  broad  belt 
f  .  irfaci  wat  ir  of  Less  than  32.  /  ><  salinity.  Tit  is  tonsists  Largely  f 

'he  remnants  of  the  diluted  coastal  belt  of  winter  which  has  beet,  carried 

i ward  by  the  same  win.:  and  turrent  systems  that  aff . h<  plumi  .  This 

less  saline  water  tends  to  obscure  the  northern  boundary  of  the  plume  at 
this  time. 

Of  particular  interest  are  the  localised  pockets  ,  f  less  saline  water 
withii  thi  river  plume  shown  in  Figures  18  and  19.  The  dif  tret  sug¬ 

gest  short-term  periodic  fluctuations  in  discharge  along  with  wind  and  tidal 
influences  which  lead  to  a  "pinching  off"  of  parcels  of  effluent  water. 

These  parcels  are  similar  to  the  'clouds'  described  by  Tully  ('°Uq)  and 
FJarli  ■  (1950).  They  are  probably  a  regular  feature  of  the  effluent  distribu¬ 
tion  in  the  sea  but  are  less  evident  during  stormy  winter  conditions  when 
mixing  is  more  ffective.  As  mentioned  above,  during  all  seasons  of  the 
year  the  ;stuary  we-  r  nters  the  nearshore  area  only  on  ebb  tid<  . 

The  plume  increased  in  areal  extent  through  late  summer  1961  and  con- 
tinued  to  spread  toward  the  southwest,  reaching  a  maximum  distance  of  about 
760  km  from  the  river  mouth  (Figures  20  and  21).  During  the  fall,  the  wind 
system  -hanged  to  southerly  and  the  low-salinity  water  was  diverted  north¬ 
ward.  As  the  river  flow  slackened  a  gradual  erosion  of  the  plume  took  place 
(Figure  22).  The  int  insification  of  the  southerly  winds  as  winter  approached 
restored  the  typical  winter  configuration  of  a  confined  low  salinity  coastal 
belt. 


Depth  Distribution  of  the  Effluent 


The  vertical  distributions  of  salinity  and  density  along  the  axes  of 
•.he  plum*  s  are  depicted  for  10  Brown  Bear  cruises  (Figure  26).  As  in  the 
areal  distribution,  the  32.5  '  /oo  isohaline  is  used  to  mark  the  Limit  of  the 
plume,  which  in  this  case  is  the  lower  boundary.  The  axis  of  the  plume, 
under  summer  conditions,  is  considered  to  be  the  line  of  Lowest  salinity 
extending  from  the  River  mouth  to  the  terminus  of  the  plume;  for  the  winter 
conditions,  the  axis  is  assumed  to  be  approximately  halfway  between  the 
coast  and  the  western  edge  of  the  belt  of  low  salinity  water  usually  defined 
by  the  32.5  L/oo  isopleth. 

The  CoLumbia  River  plume  during  the  summer  changes  the  vertical 
salinity  structure  of  the  ambient  ocean  water  as  shown  diagrammat ieally  in 
Figure  4B.  Figure  4a  illustrates  a  typical  curve  of  salinity  vs.  depth 
outside  the  plume  and  Figure  LB  shows  the  modification  of  the  salinity 
structure  brought  about  by  the  Columbia  River  effluent .  The  effluent  layer 


BROWN  BEAR  CRUISE  288(9-19  JUNE  1961) 


ACONA(OSU)  :RUISE  6)06  09-29  JUNE  1961) 


Fig.  l8--Salinity  (°/oo)  at  surface.  Brown  Bear  Cruise  288,  9-19  June 

1961  and  Acona  Cruise  6106,  19-29  June  1961  (dashed  line).  The 
distributions  for  these  two  cruises  are  shown  separated  to 
illustrate  the  effect  after  about  ten  days  of  a  strong  north¬ 
erly  wind  (l6  m  sec-l)  on  plume  movement. 


fWOWN  PFAR  CRUISI  ?‘»l  ( .’H  JUl  >  13  AUGUST  1961) 


Kig.  20--Salinity  (°/oo)  at. 
13  August  1961. 


urface,  Brown  Bear  Cruise  29 


Sigma-t  (CTf) 
Soltnity  (?••) 


Density  and  salinity  profiles  along  axis  of  Columbia  River  plume 
for  the  various  Brown  Bear  cruises,  1961-62. 
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witl  th  ;  rmanent  halt  :lin<  b<  ins  at  approximately  3<  m  iepth.  PL- 

fO-m  '  1  pth  interva.1  of  dissolved  <ygen  and  :hlor  phyll  a  maxima  Is 

Just  i  ■  1  >w  the  mixed  Layer  and  Ln  the  <  iphotic  zoni  f  r  •  h  :  n.  Fh  r< 

uv  many  fact  err  which  mu:  t  b- ■  -onsidered  when  accounting  for  th*.-  \erl  ic--l 
Li  tributi on  f  <  xyi ;en  Li  this  -egion  (P.v  kowicz,  Ln  press;  Stefans:  •  ai  1 
Richards,  in  press)  and  the  coin -idenc--  of  t he  oxygen  uni  ohlerophyl  3  maxima 
might  be  fortuitous. 

Ir.  order  to  examine  qualitatively  the  mixing  of  ;  alt  writ  r  into 
>ver  ying  freshwater,  a  (  imple  Laboratory  xj  riment  was  p  rf  rm  1.  Wind 
stress  applied  to  dyed  freshwater  overlying  salt  water  resulted  in  Internal 
waves  but  no  apparent  penetration  of  the  dye  into  th"  more  salir  wat-r. 

Wien  similar  stress  was  applied  to  clear  freshwater  overlying  dy--d  salt 
water,  ‘he  preferential  movement  of  salt  water  into  the  freshwat-  :•  .ayer 
b  evid  nt.  Figure  29  shows  diagrammatically  the  preferential  ipward 

mining.  Thin  phenomenon  has  been  observed  and  studied  in  detail  by  '  h-  rs 
(Rouse  and  Dodu  1955;  Cromwell  i960)  in  tanks  and  is  part  of  th-  ■  r.train- 
ment  process  described  by  Tully  (195^)  when  freshwater  enters  ;h<  .a.  It 
•an  be  inferred  that  the  coefficient  of  vertical  eddy  viscosity  near  th-- 
pycnocline  is  L<  ss  than  that  near  the  surface  and  thus  motions  gen<  rated 
at  the  surface  result  in  an  upper  turbulent  layer  while  there  is  r->  ia!  ively 
1 i f  t  1  e  listurbance  below  the  pycnocline  boundary.  A  parcel  of  fr-shwater 
•arried  into  the  quiet  salt  water  regime  is  buoyed  up  again  towards  the 
boundary  and  undergoes  relatively  less  mixing  than  a  parcel  of  salt  water 
■arried  into  the  upper  turbulent  regime  from  below.  In  the  ocean,  water 
from  below  is  mixed  in  the  upper  layer  of  River  effluent  and  this  mixed 
layer  expands  laterally.  Dissolved  mat  erials  in  the  river  effluent  will  be 
■arried  to  depths  primarily  by  animals,  plants  or  other  particulate  matter 
rather  that:  by  mixing  processes. 


Fresh  Water  Budget  of  the  Pi  tune 

The  balance  between  the  Columbia  River  discharge  and  the  apparent 
disappearance  of  fresh  water  from  the  plume  as  it  mixes  with  the  surrounding 
seawater  was  examined.  The  technique  employed  is  based  on  the  assumption 
that  the  bae  kg round  ocean  water  has  a  salinity  of  32.5  l/°°  and  that  sali¬ 
nities  in  a  water  column  less  than  32-5  /oo  result  solely  from  a  mixture  of 
riv  -r  water  and  the  ba  •  kg round  water.  If  this  fresh  water  "fraction”  in  the 
wat  -r  -o.lumn  could  be  separated  from  the  ocean  water,  its  equivalent  height 
is  given  by: 

A 

„  _  32.  d  Z*  —  o.)  S  dZ 


wh-'t~"  H  ir  th-'  height  in  meters  of  fresh  water,  S  the  observed  salinity. 

32.3  th  salinity  of  the  ocean  background  water,  Z  the  depth  and  Z*  the 
d-'pth  of  th-'  32.5  <’/oo  isohaline.  The  defects  in  this  appro&-c4aare  (a)  the 
a  Li  :  of  ‘h-’  background  ocean  water  varies  a  little  from  32.5  /oo,  and 
(b)  cth-'r  possible  sources  of  dilution,  e.g.,  precipitation,  are  ignored; 
how'v->r.  it  is  felt  that  the  technique  described  above  is  useful  as  a  first 
'Approximation  in  estimating  the  volume  of  river  water  in  the  plume.  Ketchum 
and  Keen  (l  v'5)  and  Tully  and  Barber  (i960)  have  used  similar  approaches  ;i 


«. 
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tudl  f  listrlbutioi  f  *  and  rain  wat  r  I  ar  areas  1  t) 

'  •  -sin. 

T‘  >  conputat  ior.i  d  .  Tib  a  abov-  wen1  ••on.pb't'  i  is:  ;h-  witut 

. ■  ry  '  ■  ■  :  :  '.x i  t  h  tl  [  i  for  :ru3  Lsi  1  -  Pal 

Larg"  -si  1  plots  of  the  vertical  i ribution  of  salinity  v»'n  trad-  • 
each  st at  ion  a: id  the  value  of 

Z* 

^ j .  5  z*  -  J'  saz 

was  btai  i  is  ing  a  planiraeter.  Phis  .  wa  '  1  n  livid  d  by  •.  .  -  to 
btaii  th  qui valent  fr  shwat  r  1  '  ht  Ln  n  t  rs.  N  un  irical  Lnl  ratios 

was  als  sarrj  it  for  mos  t  of  th  sruises.  Hie  different  It  heights 
betw  ;  tw  method  was  ss  s  1  ha  j  rcent  and  th  sompul  ■  ;  r 
for  del  rn  Inin)  quival  nt  fits  1  wat  sr  heights  Ls  satisfact  ry . 

:■  : harts  >f  quivalenl  fresh  water  heights  wer  prepared  for 
.0  cruis  •  .  and  th  -  sippsirent  freshwater  volume  was  determined  by  planimetry 
n  ach  f  these  iharts.  Ph<  results  f  thess  ;omp  itat  ions  ar  showi 
rat  •  nd  :ontoured  iharts,  which  typi fy  winter  and  s  unm  ■  ;oi  lit  Lons  . 

ar  ••  nt  1  for  four  sruises  (Figures  -  to  33) •  From  January  -t  to 

Jai  u  ry  .  th  av  rags  lai  y  lisappearance  or  "loss"  of  fr  shwat  ;r 
(i.  mixing  shwat  r  witl  .  awat  r  intil  it  ig  n<  Long  r  deti  itabl 

by  the  32.  >  :rit  rioi  )  was  fot^t  i  to  be  6.0  x  10^m~aay~  ,  a  n  an  annual 

f]  in  water  f  6»  CB  sec  ’  (232,000  ft  s  If  ::  Ls  assumed 

hat  tl  r  Ls  no  net  iccumulatioi  f  freshwater  or  increas  salt  :ontent 
..  t  m  .  irrounding  ths  plum  .  tl  ■  the  mean  annual  river  dis  :har 
1  e  u  "  oss"  of  freshwat  r  from  the  plum  .  rh  p  um  "loss" 

6,600  agr  wel  wit  in  annua  lischar  of  7, 300  ursec  . 

agr  :m  nt,  although  perhaps  fortuitous,  supports  the  sel  jtion  of  th 

32.  oo  I  oha  Lns  tl  boundary  of  th  j  un  ■ .  rh  valu  s  showi  Ln  Pabl 

U  ar-’  from  calculations  made  reg i > ■ ‘t ing  pr "cipitation  and  -vaporntion  and 
• :  :on1  -ibution  of  soastal  str  ami  .  e  negative  value  Ln  the  "  oss"  of 

fn  hwat  ■  vs  Lui  Ls  probably  a  r  suit  f  the  :ombined  sffects  of  the  rapidly 
lees  Le rat ing  river  runoff  Ln  mid-June  (Fig.  2)  and  the  wind  change  showi 
Figure  lb.  The  excess  apparent  freshwater  volume  during  this  period  might 
also  be  attributed  to  the  return  southward  of  the  Columbia  ar.d  coastal  river 
effluent  which  were  previously  discharged  and  lying  to  the  north  cl-  the 
Columbia  Ri\  r  mouth.  If  the  average  daily  "loss"  of  squivalent  fr  hwat  • 
volume  is  normalized  tq  th"  area  cf  the  plume,  it  is  found  to  be  npproxi- 

bm  tel;  ;m/day  m  .  The  volumes  of  freshwater  obss  rved  f< . ich 

p-riod  of  the  survey  and  those  for  the  daily  river  discharge  are  presented 
in  Figure  tU. 


VII.  MAJOR  PHYSICAL  PHENOMENA  AFFECTING  THE  DISPERSION 
OF  THE  EFFLUENT  IN  THE  SEA 

General  Description 

Aft  -r  initial  mixing  in  the  estuary  cf  about  1  part  river  water  with 
•  par  .  th  ffluenl  nters  the  sea  as  a  "Jet"  of  3 ow  salinity 


Fig.  32 --Equivalent  height  of  fresh  water  in  meters.  Brown  Bear  Cruise 
293,  1*+  September-20  October  1961. 


wit  r  where  it  continuin'  to  bo  mixed  by  hydraulic,  tidal,  and  surf-act  lor: 
processes  if.  t If  nearfhot’"  ntva  --xtending  '.>0  km  from  the  River  mouth.  It. 
this?  area  lateral  .and  vertical  mixing  processes  .•rent**  :i  ’.in*  :  our  q- 
relatively  freshwat  *r,  a  i  Lxture  of  about  pari  river  wal  *r  o  parts 
;  ■‘••iw.at  ■■!•  about  40  km  In  length  lying  >0  to  >0  km  seaward  of  tit*  mouth.  A. 
tit**  '. ow  salinity  mixture  of  river  and  salt  water  accumulates  loyally  on  top 
of  denser,  mor*  saline  ocean  wa.tr,  a  surface  pressure  gradient  is  main¬ 
tained  which  <  aids  to  move  the  freshwat  ;r  away  from  t  h  sourc  .  P 
partially  mixed  effluent  niters  into  the  oceanic  environment  where  different 
processes-  of  mixing  and  diffusion  predominate.  In  the  offshor-  region  ad¬ 
vert  ion  by  wind- induc'd  and  geostrophic  currents,  and  eddy  diffusion  play 
major  roles-  in  dispersing  and  eroding  the  plume  illustrated  diugrammat  i  cn  1  ly 
for  s ummer  conditions-  in  Figure 


Curas  at : 


The  currents,  deduced  from  t-he  geopotential  topography  follow  a 
seasonal  pattern;  they  set  south  both  inshore  and  offshore  at  b  to  •>  miles 
per  day  in  the  summer  and  north  at  about  j  miles  per  day  along  t  tv.*  const 
during  the  winter.  From  the  month- to-month  details-  of  -urr-nt  direction 
atui  speed  it  is-  possible  to  determine  the  average  direction  and  magnitude 
of  plum**  movement  caused  by  these  uirrents.  In  some  areas  of  the  region 
large  -'ddies  frequently  complicate  the  pattern  and  during  most  periods-  only 
a  rough  average  can  be  estimated.  For  example,  the  southerly  setting  cur¬ 
rent  off  the  coast  of  Oregon  in  the  summer  of  1961  was-  about  h  cm  sec-';  the 
range  varying  locally  from  l  to  00  cm  sec- ' .  Direct  ion  is-  also  variable 
(of.  Figures  ‘*  to  lit).  The  apparent  eddies  might  be  associated  with  internal 
waves-  as  mentioned  earlier.  Geostrophic  currents  are  computed  on  the  as¬ 
sumption  of  a  steady  state  distribut  ion  of  mass-;  whereas  transitory  •  * f f oot  s 
such  as  internal  waves  associated  with  tides  or  tin*  passage  of  local  storms: 
or  intrusion  of  river  effluent  may  cause  t-he  local  mass  structure  to  depart 
s  igni  f  icantl.v  from  steady  state  conditions  and  lead  to  d  isorepanc  ies-  in 
computed  currents-.  Recog.nir.ing  t hose  limitations,  the  currents  have  beer, 
computed  for  various  periods  of  the  year  and  the  results  used  along  with 
wind  un-rents  to  estimate  the  resultant  movement  of  tin*  Columbia  River  water. 


As  an  initial  approximation  the  geostrophic  current  has  been 
set  l80°T  at  b  cm  sec-'  it:  the  summer  and  000°T  tit  6  cm  sec- 


er-t- limited  to 


in  t  tie  winter. 


Wind  Transport  of 


As  mentioned  earlier,  the  seasonal  wind  patterns  of  the  area  off  the 
•oasts  of  Washington,  Oregon  and  Northern  California  consist  of  predominant ly 
northerly  winds  of  .about  6  m  sec-1  in  summer,  and  southerly  winds,  frequently 
of  gale  force,  in  winter.  The  wind  drives  and  mixes-  the  freshwater  dis- 
•harg'd  into  tin*  open  ocean,  greatly  affecting  t-he  distribution  of  the 
Co  lumhi.u  R  i  vet-  off  luent  . 


The  Internet  Ion  between  the  wind  and  surface  waters  has  been  the  sub- 
.]■'  -t  of  cons  Iderable  controversy.  A  thorough  understanding  of  this  inter¬ 
act  lor  lias  yet  to  emerge  from  t-he  numerous  studies.  The  velocity  of  the 
urfa  *e  drift  current  is  a  fund  ion  of  the  applied  surface  wind  st  ress  and 


the  eddy  viscosity  the  water.  Becaus<  of  the  relatively  small  -idy 
'•is  .*»':•  i t  y  rear  the  shallow  sharp  pycnocline  in  the  plume  area,  the  momentum 
I:.”:-  will  he  confined  to  a  shallower  layer  in  the  plume  area  than  outside 
the  pi  time  in  the  more  homogeneous  water.  The  program  for  •omputing  the  wind 
■urr  at  is  die  sussed  below. 

The  driving  stress  of  the  wind  on  the  s."i  surface  -ausee  a  surface 
•urrent  which  is  directed,  approximately  45  degrees  to  the  right  of  the  wind 
Li  th  northeri  hemisphere  (Ekman  L902;  Bossby  and  Montgomery  193  )•  The 
water  speed  '  s-  related  directly  to  the  wind  speed  by  the  empirical  relation 
derived  by  Ekman: 


V  = 


0.0127 

/s  i  n  Q 


(1) 


wher  W  is  the  wind  speed  and  Q  the  latitude.  The  wind  factor.  0.0127, 
has  been  examined  by  a  number  of  investigators  and  a  recent  analysis  of 
over  r7,000  observations  supports  this  value  (Budinger,  unpublished).  The 
•quation  is  applicable  only  to  the  surface  wind  current  and  tells  little 
about  the  transport  by  wind-induced  currents  across  a  vertical  surface  of 
unit  width.  Tills  trans-port  (.'ailed  the  Ekman  transport)  is  related  to  the 
wind  stresi  by  Int  ‘gration  of  the  equations  of  motion.  In  open  water  the 
components  are 


Ex  Ty/f 

and  (2) 

E.V  -  T  x/f 

where  E\  and  Ey  are  the  conn l  and  meridional  Ekman  mass  transports;  tx 
and  Ty,  the  surface  stresses;  and  f,  the  Coriolis  parameter.  Tlte  net 
transport  is  independent  of  any  change  in  eddy  viscosity  with  depth  and 
the  mean  transport  is  directed  (in  the  northern  hemisphere)  9^  degrees  to 
the  right  of  the  wind  direotion. 

Th-3  surface  stresses  can  be  calculated  using  the  theory  that  the 
stress  is  proport  tonal  to  the  square  of  the  wind  velocity: 

!  T  I  =  CD  |  w|-\  (3) 

where  p  is  the  air  density,  W  the  wind  velocity  at  10  m.  CL.  is  a  dimen- 

^  d  *  19 

sionless  drug  coefficient,  and  W  is  the  wind  velocity  at  10  m.  The 
validity  of  this  theory  and  the  magnitude  of  the  drag  coefficient  have 
been  under  investigation  for  many  years .  Recently,  Wilson  (i960)  examined 
the  results  of  nearly  50  independent  investigations  and  evaluations  of  the 
stress  law,  and  concluded  that  the  mean  value  of  the  drag  coefficient  is 
0.0024  for  strong  winds  and  0.0015  for  light  winds.  The  two  different  co¬ 
efficients  have  been  attributed  to  an  increase  in  the  roughness  of  the  sea 
surface  as  the  wind  increases.  The  agreement  among  the  authorities  quoted 
was  fairly  satisfactory  for  strong  winds,  but  less  so  for  light  winds. 

Tiie re  is  no  universal  agreement  on  the  velocity  or  range  of  velocities 
demarking  strong  from  light  winds  in  computing  drag  coefficients.  Deacon 
and  W-'bb  (bH'2)  summarised  careful  observations  made  under  conditions  of 
neutral  air  stability  and  found  a  range  of  drag  coefficient  values  below 
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0.00.'.  They  suggest  that  th<*  drag  coeffi  eient  probably  increasm 
with  wind  speed .  This  relatively  .  Light  increase  is  described  by 

0  (l.OO  +  0.07W)  x  1C'-'’ 

where  W  :  the  wind  speed  in  m  sec  1  at  10  meters.  Thus,  the  drtig  co¬ 
efficient  for  tiie  average  winds  of  8  tti  see-1  during  the  summer  would  be 
0.0016  and  for  winds  of  20  bus  so” 1 ,  0.0024.  Sheppard's  ( )  equation 
gives.  0.0010  and  0.0024,  respectively. 

The  men::  stress-  exerted  by  a  fluctuating  wind  depends  not  only  on 
the  drag  :oeff icient  but  also  on  the  value  to  be  used  for  W  .  Practi  tally, 
winds  must  be  taken  from  charts  and  hence  have  been  smoothed  over  con¬ 
siderable  periods  of  time  (monthly  or  daily),  and  as  (W)'  is  always  small  s 
than  (W-).  a  larger  drag  coefficient  is  necessary  to  give  the  same  stress 
value  than  would  be  the  case  if  winds  were  measured  continuously.  Further¬ 
more,  the  values  of  the  drag  coefficients  estimated  from  different  time- 
averaged  observations  will  vary,  and  at  present  there  is  no  means  by  which 
wind  can  be  accurately  translated  into  water  movement  under  all  conditions. 
From  a  review  of  the  work  of  Taylor  (.1016),  Ekrnan  (1928),  Rossby  and 
Montgomery  (’.98^),  Palme n  and  Luurila  (1938),  and  Sverdrup  and  Fleming  (194  i) 
and  tin'  results  of  the  use  of  different  coefficients  under  various  natural 
conditions,  a  value  of  0.0024  for  his  drag  coefficient  appears  to  be  the 
most  reasonable  for  use  in  the  Columbia  River  effluent  area. 

The  wind  speeds  in  t'  .  program  were  computed  from  twice  daily  atmos¬ 
pheric  charts.  The  effect  of  atmospheric  stability  will  be  such  as  to 
reduce  the  drag  coefficient  for  stable  conditions  and  increase  the  coefficient 
for  unstable  conditions;  however,  the  exact  relationship  has  not  been  deter¬ 
mined  . 

If  the  wind  field  is  known,  its  effect  on  the  water  can  be  comput'd 
using  the  -quations  (l)  and  (_')•  It  is  possible  to  derive  the  surface  winds 
and  thus  the  surface  stresses  from  t  lie  geos  trophic  winds,  which  in  turn  are 
a  function  of  the  atmospheric  pressure  field:  t  a  VT~  a  A  F  .  The 
pressures  are  available  from  the  meteorological  reports.  This  method  ,9’ 
determining  the  surface  stresses  is  similar  to  that  used  by  Montgomery  ( 193!' ) 
and  Fofonotf  (l9o0).  Fofonoff,  applying  the  theory  that  total  steady- state 
transport  of  mass  in  the  ocean  depends  primarily  on  the  curl  of  the  wind 
stress  noting  on  the  surface  of  the  ocean  (Sverdrup  1947;  Stommol  ip48;  Munk 
19a0),  computed  the  transports  in  the  North  Atlantic  and  North  Pacific  Oceans. 
He  used  monthly  means  of  sea- level  pressure  (averaged  from  twice-daily 
observations )  and  plotted  these  on  a  10-degree  diamond  grid.  The  atmospheric 
pressures  were  taken  from  the  observations  made  by  the  Extended  Forecast 
Section  of  the  U.S.  Weather  Bureau.  It  was  found  that  these  data  were  in¬ 
adequate  for  the  Columbia  River  area  because  day-to-day  variations  as  well 
as  small-scale  space  variations  arc  important  in  this  problem.  I’e  obtain 
adequate  information  it  was  necessary  to  determine  values  of  a :  osrher c 
for  si  hi  rid  points,  shown  Ln  Figure  36,  from  00CX  a  d 
Or  •  .  ■■■  harti  of  atmospheric  pressure  prepare >a  bj  tl  Loct  w  '  tr  ir  au 
•it  the  S'  at  t  i. '-Tacoma  Inf-mat  .Iona  L  Airport. 
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Fig.  36--Monthly  averages  of  Ekman  transport,  1961,  converted  to  velocity 
(cm  sec"T)  of  the  water  column  0  to  to  meters. 


Pressures  were  recorded  for  each  of  the  eight  grid  points  for  the 
period  of  January  1961  through  iy62  to  present.  The  University  of  Washington 
IBM  709  computer  was  programed  to  compute  the  meridional,  zonal  and  result¬ 
ant  surface  wind,  surface  drift  current,  wind  stress,  and  Ekman  transport 
from  the  twice  daily  estimates  of  the  pressure  field  off  the  Washington  and 
Oregon  coasts.  The  computation  scheme  entailed  determining  the  geoetrophic 
wind  velocity  aloft.  This  wind  vector  was  then  rotated  15  degrees  to  the 
left  (Montgomery  1935)  of  its  downwind  direction  and  reduced  by  30  percent 
(Gordon  1950).  The  result  of  this  operation  is  the  10-meter  wind  (the  W 
applicable  to  the  equations).  The  meridional  and  zonal  components  of  this 
wind  were  then  used  to  determine  the  surface  current  using  equation  (l). 

The  direction  of  this  current  was  assumed  to  be  45°  to  the  right  of  the 
surface  wind  vector. 

The  surface  stress,  proportional  to  the  square  of  the  wind  speed,  was 
also  calculated  using  the  stress  relationships: 

T  x  =  p  C_  U  h2  +  V2 
'a  D 

ry  =  pa  CD  V  M  +  V2  , 

where  U  and  V  are  the  zonal  and  meridional  components  of  W,  respectively. 

The  air  density,  p  ,  was  taken  as  1.22  x  10"  3  gm  cm"3,  and  the  drag  co- 
efficients  Cp,  as  2.4  x  10“3.  if  the  wind  speeds  are  in  units  of  cm  sec”  , 
the  stress  is  in  dynes  cm"^.  The  Ekman  transport  was  computed  using  the 
computed  surface  stresses  and  the  Coriolis  parameter  according  to  (2). 

Because  the  eddy  viscosity  becomes  very  small  at  the  lower  boundary 
of  the  plume,  the  depth  of  frictional  resistance  of  the  Ekman  layer  is 
assumed  to  coincide  with  the  bottom  of  the  plume,  at  about  40  m.  Therefore 
the  movement  of  a  water  particle  within  the  plume  due  to  the  wind-induced 
Ekman  transport  can  be  calculated  from: 

I  -  I  |  -  I 

.  ,  '  V  ‘  [  '  °'57 

where  '  •  is  the  horizontal  displacement  of  a  1  m  x  40  m  column  in 

nautical  miles  per  12-hour  period  and  j  E  is  the  Ekman  transport  in  metric 
tons  rrT^sec-^.  The  resultant  monthly  averages  of  Ekman  transports  are  pre¬ 
sented  in  Figure  3 6  for  the  year  1961.  These  vectors  clearly  show  the 
seasonal  variation  of  the  wind  effect  on  the  surface  transport.  The  net 
Ekman  transport  for  1961  is  south.  Power  spectrum  analyses  and  further 
studies  are  planned  for  the  values  of  Ekman  transport. 

The  period  before  a  steady  transport  is  reached  for  a  particular 
wind  stress  has  been  estimated  to  exceed  5  pendulum  hours  and  is  probably 
in  the  order  of  one  day  (Ekman  1905;  Von  Arx  1961).  By  integrating  the 
transports  during  each  cruise  or  between  mid-times  of  successive  cruises, 
the  magnitude  and  direction  of  the  wind-induced  movement  can  be  obtained. 

Then  the  vector  sum  of  the  geostrophic  current  and  the  Ekman  tra.nsport 
current  gives  the  translation  of  the  plume  for  the  period.  Figure  37  shows 
the  predicted  axes  for  Cruises  288  (June),  290  (July),  and  293  (September- 
October  1961).  A  constant  4  mile  day"^  current  setting  south  was  used, 


Fig-  37--Predicted  axis  and  outward  extent  of  plume:  (A)  Brown  Bear 
Cruise  288,  9-19  June  1961,  and  Acona  Cruise  6106,  19-29  June 
1961  (dashed  line);  (B)  Brown  Bear  Cruise  290,  6-25  July  1961 
(c)  Brown  Bear  Cruise  293 >  lV  September-20  October  1961. 


:xeept  r  a  j  ri  d  pr  i<  Ling  Jrui:  <  )(  when  ther  was  nc  :urr  nl  in  at 
area  of  di\  rgenc  (Pigur  ) ) .  During  Jum  whei  the  lisehargi  i:  the 
reatest  ai  Ekma  Lay  sr  f  .  m<  t  rs  thd  :kn<  ss  was  assumed;  for  (  ther 

periods  •  i:  ■  layt-r  war  us:  .L’Ti"d  to  t  •  40  I'.-'t-rs, 

Pd  ir  A  sh  ws  ■■  rved  aiovem  nt  of  the  plum  luring  a  10-day 

period  In  June.  The  origin  for  the  initial  v<  'tor  was  4(  to-:  si  of  the 

Rivi  r  mi  ith.  A  strong  n  rtl  r  y  storm  (K  m  .  ■  ■  • )  :ommi  •  :ed  luring  the 

la:  r  part  of  Cruis  •  2! (9-19  June)  and  abated  during  the  firs t  part  V 
ACONA  Irui  6l(  (  -  >-!  1  J  a  • N .  The  predicted  :enter  f  th«  tell  is  the  re¬ 
sultant  f  •  h  \  :tor  Iditi  >n  of  th  ostrophic  current  and  wii  L  driven 

:urren1  f  r  the  ;  riod  L0-2C  June.  The  projected  /oo  isohalin  was 
onstrucl  • i  ..  th  sam  .  :tor.  Figure  37B  shows  the  predicted  axis  of 
:  r.  •  plum-  oar  -  d  on  the  geos  trophic  and  wind  currents  from  Cruise  2Ht  (June) 
to  Cruis  (July).  Pi;  predi  :1  1  axis  for  September-October  1  -  1  (Pig. 

37c)  was  prepan  i  .  aying  out  th  r  sultant  vectors  for  the  middle  and 
Lower  pridt  f  Figure  with  an  origin  ho  km  off  the  River  and  on  the  axis 
f  Cruis  290,  r  spectively.  Somewhat  better  agreement  between  the  actual 
and  predicted  axes  for  the  southern  section  of  the  plume  would  be  obtained 
if  a  grid  for  Ekman  transport  had  been  :onstructed  farther  south — the  wind 
current  rear  40°N  was  northwest  during  the  latter  part  of  the  prediction 
rlod  fr  -:\  .  •  July  to  2  Octol  L96I.  The  predicted  axes  and  water  move- 
are  in  fair  agreement  with  the  observed.  It  should  be  str  ssed, 
however,  thal  axes  and  the  cell  trajectories  are  subject  to  errors  in 
:hoic  if  computation  par am  ters  (i.e.,  layer  depth,  drag  coefficient, 
av  raging  period,  Lir  :tion  and  speed  of  the  geostrophic  currents). 


VIII.  DIFFUSION  OF  RIVER  WATER  IN  THE  SEA 


J-  n-'ral  R  ‘marks 


Although  the  vector  sum  of  the  geostrophic  current  and  the  velocity 
d-. -rived  from  Ekmar.  transport  computations  roughly  describes  the  advert  ion 
of  -f fluent  water,  the  model  for  analysis  and  prediction  is  incomplete  un- 
Les!  diffusion  is  also  considered.  The  yearly  average  rate  of  detectable 
effluent  disappearance  (mixed  beyond  recognition  by  the  measuring  technique- 
approximate  iy  part  f  River  water  to  1,000  parts  seawater)  was  computed 
as  about  <  .0  x  LCr-  m3  day  ”  .  This  figure  agrees  veil  with  the  coapui  1 
a\  >rage  annual  River  discharge.  It  appears  that  the  effluent  Wat  -r.  Lying 
on  the  surface  become  mixed  vertically  with  seawater  from  below  the  upper 
ha  Iodine  by  a  preferential  upward  movement  of  deep  water,  and  laterally, 

■is  the  plume  spreads,  with  background  water  of  32.5  ''/oo.  Tills,  results  in 
a  slight  increase  in  thickness  (about  10  m  in  280  Kilometers.)  and  width  of 
the  pan  cake - shaped  effluent  with  time  and  distance  from  the  River  mouth. 

Vertical  and  lateral  diffusion  processes  in  the  sea  have  proved 
difficult  to  describe  and  analyze,  and  previous  attempts  to  describe  the 
diffusion  of  an  effluent  in  the  sea  are  not  wholly  satisfactory  for  this 
syst  an  (Takano  1954,  L954a,  1955}  Bamada  1059;  Gifford  1959;  Brooks  L960, 
and  Defant  bVl).  For  most  analyses  it  is  necessary  to  determine  the  hori¬ 
zontal  and  vertical  coefficients  of  eddy  diffusion,  K  .  K  .  and  K  .  These 

x  y  z 
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■oefficients  are  used  in  a  diffusion  law  of  one  form  or  another  (Pearson 
i960;  Waldichuk  1963).  Reviews  of  these  coefficients  have  been  given  by 
Okubo  (1962)  and  Ichiye  (1962). 

As  mentioned  previously  the  problem  of  discharge  into  the  ocean  might 
best  be  separated  into  an  examination  of  the  nearshore  tidal  and  hydraulic 
regime  (less  than  50  km  from  the  estuary)  and  the  offshore  wind  and  current 
regime.  Work  has  just  begun  in  the  more  complicated  nearshore  environment, 
and  no  results  are  as  yet  available. 


Vertical  Mixing 


The  river  effluent  is  confined  to  the  first  few  tens  of  meters  (Fig. 
4  and  29)  and  is  mixed  with  the  ambient  seawater  by  upward  transport  of 
properties  from  the  deeper  water  and  horizontal  turbulence.  Three  methods 
were  used  in  an  attempt  to  evaluate  the  vertical  coefficient  of  eddy 
diffusion. 


METHOD  I--The  vertical  salinity  section  along  the  plume  axis  for 

AS  ^2c 

Cruise  288  (June  1961;  Fig.  26)  was  used  to  compute  K  from  u  —  =  K  — g. 

.  Z  OX  Z  C3Z 

The  ratio  Kz/u  was  determined  at  10-meter  depth  increments  down  to  40  meters 
(the  approximate  bottom  of  the  plume)  for  6  locations  30  km  apart  along  the 
axis  starting  50  km  from  the  source.  The  values  for  K  /u  were  not  consist¬ 
ent,  often  changing  in  sign  and  magnitude.  An  average  of  the  positive  values 
of  K„/u  was  approximately  0.3  cm.  Assuming  a  realistic  value  of  10  cm/sec 
for  u,  average  is  3  cm^  sec“^.  A  smoothed  composite  vertical  salinity 
distribulion  approximating  a  combination  of  the  profiles  from  Cruises  288 
(June  1961),  290  (July  1961),  293  (September-October  1961),  and  30C  (June 
1962)  was  then  used  to  compute  K  /u.  It  was  hoped  that  the  smoothing  would 
remove  some  of  the  discrepancies  in  the  values  of  K  /u;  however,  the  mean 


value  of  K  of 
z 


2  -1 
cm  sec 


was  averaged  from  a  wide  scatter  of  values. 


METHOD  II- -The  second  method  was  based  on  the  Fickian  equation  of 
diffusion 


dM  =  -K_  dfc] 

It  Z  dn 

where  the  mass,  M,  of  a  substance  diffusing  through  a  unit  area  in  time, 
t,  is  proportional  to  the  concentration  gradient  normal  to  the  unit  area. 
The  vertical  eddy  coefficient,  K  ,  can  be  evaluated  from  this  formula  if 
estimates  can  be  made  of  the  change  of  mass  with  time  within  one  water 
mass.  During  Cruise  309  (September  1961),  a  series  of  shallow  stations 
(4TA  to  47H)  was  made  in  a  manner  which  allowed  evaluation  of  the  change 
in  mass  of  the  water  column  above  the  halocline  over  a  period  of  time.  A 
24-foot  aviator's  parachute  was  used  as  a  sea  anchor  and  the  ship  occupied 
hydrographic  stations  at  approximately  3-hour  intervals  for  a  total  of  21 
hours  during  which  the  ship  drifted  south  at  15-25  cm  sec-1.  Figure  38 
shows  the  vertical  distribution  of  salinity  at  stations  47A  to  47H.  This 
figure  shows: 
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Fig.  38--Verticnl  distributions  of  salinity  at  Brown  Bear  drifting  stn 
tlon  1+7  during  Cruise  109»  20  June-9  July  19^2. 
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METHOD  III--When  the  plume  is  still  growing. its  volume  increases 
downstream  and  the  change  of  volume  from  one  section  to  the  next  may  be 
converted  to  an  addition  of  water  of  salinity  j2,5  l/oo  per  unit  time.  An 
apparent  upward  speed  can  be  computed  by  dividing  this  volume  addition  per 
unit  time  by  the  average  surface  area  of  the  two  sections  used  to  determine 
the  change  in  volume.  The  plume  near  its  origin  during  Cruise  288  (Fig. 

18)  was  divided  into  three  adjacent  solid  rectangles  oriented  normal  to  the 
axis.  The  bottom  of  these  adjacent  sections  was  the  32.5  /oo  surface. 
Assuming  that  all  of  the  salt  was  added  to  each  section  by  advection  from 
below,  the  mean  salinity  for  the  next  section  can  be  calculated  and  com¬ 
pared  with  the  integrated  salinity  from  the  observed  values.  The  results 
of  this  technique  and  the  vertical  speeds  are  shown  in  Table  6.  The 
vertical  upward  speeds  decrease  with  an  increase  in  distance  from  the  mouth 
of  the  t'iv-‘r.  The  vertical  speed  of  about  1.5  meters  day-'  (l.7  x  10--^  cm 
sec-  1 )  is  in  the  order  of  magnitude  of  speeds  estimated  for  upwelling  in 
'th"t-  situations  by  McEwen  ( 1 93^+ ) ?  Sverdrup  and  Fleming  (194 1),  Saito  (1951) 
and  by  Hidaka's  (195*0  quantitative  conclusions  based  on  a  theoretical  model 

To  compare  the  results  of  this  method  with  those  from  METHODS  I  and 
II  which  give  the  apparent  change  as  due  to  diffusion  (i.e.,  in  terms  of 

an  eddy  ’oeff ieient ),  the  term  for  vertical  advection  ,  ^  _dS  \  in  the  distri 

button  >f  variabl  s  equation  was  equated  to  '  bZ  ' 

.  Assume  K.  is  constant  and  W  from  METHOD  III  is  1.7  x  10  J  cm  sec 

^  * 1  *J  dS  ^2  S  -2 

Typical  values  for  — -  and  -rrp-  for  the  area  of  observation  are  3  x  10 
7  .  dZ 

l/oo  -m  and  4  x  10“ 1  L/oo  cm-1  ,  respectively.  The  coefficient  K  computed 
by  equating  the  advection  term  to  the  diffusion  term  is  1.4  cm‘°  sec-*. 
METHODS  I  and  II  give  coefficients  of  eddy  diffusion  between  2  and  11  cm' 
sec-*.  The  evaluation  of  the  change  of  upward  velocity  with  distance  from 
the  origin,  horizontal  advection,  and  stability  is  being  pursued. 

OTHER  METHODS — The  closely  spaced  observations  shown  in  Figure  26 
and  more  recent,  sections  through  the  plume  edge  provide  a  basis  for  the 
estimation  of  Kz  and  its  relation  to  stability  and  distance  from  the  river 
mouth.  The  undulations  in  the  halocline  must  be  smoothed  before  computa¬ 
tions  based  on  the  equation 

u  2§  =  Kz 

dX  az7 

will  be  valid.  These  computations  have  not  been  made  as  yet  . 

Measurements  of  horizontal  and  vertical  diffusion  of  planted  dye 
and  radioactive  materials  are  planned  and  estimation  of  the  coefficients 
will  be  made  using  the  techniques  outlined  by  Sutton  (l953)>  Crank  (1958), 
Pritchard  and  Carpenter  (1980),  and  Okubo/Mtz), 


Horizontal  Diffusion 


Horizontal  coefficients  of  eddy  diffusion,  K  ,  and  K  ,  have  been  shown 
to  vary  with  a  scale  of  the  phenomenon.  This  scale  is  sometimes  considered 
t.o  be  the  distance  between  neighboring  particles  or  the  size  of  an  eddy;  or 


TABLE  6 


COMPUTATION  OF  VERTICAL  VELOCITIES  FROM  CONTINUITY 


S2V2  =  S^V,  +  (V2  -  V,)(32.5)  (1) 


Subscripts  refer  ction  numbers. 

S  =  mean  salii.  ty  (  / 

V  =  volume  (cu.  in. ) 

A  =  area  (sq.  m. 

W  =  vertical  velocity  (m/day) 

U  =  mean  current,  taken  to  be  5  nm/day 

X  =  translation  distance  due  to  U, 


10 

run.  (width  of  se 

ctions ) 

Mean  Salinities 

Integrated 

Computed  from 

Computed  from 

W-1 

Region  No. 

from  contours 

equation  (l) 

equation  (3) 

m  da 

1 

31-57 

1.72 

2 

31.69 

31.68 

1.47 

3 

31.70 

31.77 

32.31* 

0.162 

4 

31.65 

31.71 

34.67* 

*  The  W  calculated  from  two  preceding  volume  segments  is  used  to  pre¬ 
dict  the  S  of  the  next  section;  since  W  decreases  with  distance  from 
the  source  of  the  freshwater,  this  W  is  too  large  for_the  situation 
existing  in  the  next  region,  and  hence  the  resulting  S  is  too  large. 
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the  size  of  the  water  body  in  more  restricted  waters  (see,  e.g.,  Sverdrup 
19^6;  Stommel  194-9)-  After  the  scale  of  phenomenon  is  defined  various  "laws" 
have  been  used  to  relate  the  eddy  diffusion  coefficients  to  the  length 
parameter.  The  relation  usually  takes  the  form  K  a  LX  where  x  lies  between 
1  and  4/3  depending  on  the  assumptions  and  arguments  used.  An  extensive 
examination  of  the  literature  for  values  of  K  calculated  from  distribu¬ 
tions  of  properties  in  nature  (Coachman  and  B^dinger,  unpublished)  revealed 
that  a  value  of  x  between  1  and  1.2  best  fitted  the  observations,  though 
there  is  always  considerable  variation  in  this  type  of  calculation.  It  is 
obvious  that  the  coefficient  will  be  subject  to  the  choice  of  the  scale  of 
the  phenomenon  and  in  many  cases  this  choice  is  quite  arbitrary. 

Two  techniques  were  used  to  calculate  K  and  K  : 

x  y 

METHOD  I--Over  most  of  the  plume  beyond  40  to  50  km  from  the  river 
mouth  it  can  be  assumed  that  advection  along  the  axis  is  balanced  by  lateral 
and  vertical  diffusion: 


Kz 


d^C 

dZ- 


where  U  is  the  velocity  along  x,  the  axis  of  the  plume;  C  is  the  concentra¬ 
tion  of  freshwater,  and  K  and  K  are  constant  eddy  coefficients.  Approxi¬ 
mate  solutions  for  this  equation  have  been  given  by  Roberts  (1923).  The 
general  solution  for  a  continuous  point  source  introduction  of  material 
which  is  advected  in  x  and  dissipated  through  non-isotropic  mixing  is: 


X 


Q 


+  yf  + 


4tt  ■(  K  K  K 

L  x  y  z  J  L  Kx  Ky 


: 


P 

i 

1 r 


If  diffusion  in  the  X  direction  is  small  compared  to  the  advection;  that 
is,  if 

l  Ky _ Kz_  J 

jL 

Kx 


is  negligible,  the  solution  reduces  to 


C 


Q 


4rr  (Ky  Kz)2  x 


U  /  yf  +  zf\ 
¥x  VKy  Kz  J 


J 


Values  for  Ky/Kz  can  be  obtained  from  this  equation  by  applying  data  from 
various  cruises  and  the  river  discharge.  In  order  to  apply  the  equation, 
certain  assumptions  regarding  the  quantities  U,  C,  and  Q  were  necessary. 
First,  the  maximum  concentration  of  freshwater  on  a  chart  similar  to 
Figures  30  to  33  was  taken  as  the  axis  of  the  plume,  the  x-axis.  The  value 
of  U  was  estimated  from  the  distance  along  the  axis  to  the  edge  of  the  plume 
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i 
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divided  by  the  number  of  weeks  of  accumulated  water  in  the  plume.  Thi. 
method  leads  to  advection  values  which  are  somewhat  lower  than  those  deter¬ 
mined  from  the  wind  and  geostrophic  currents;  however,  the  technique  used 
to  evaluate  Ky/Kz  is  not  sensitive  to  the  value  of  U.  The  concentration  of 
freshwater  in  the  plume  (C)  was  defined  as  the  dimensionless  ratio 
Equivalent  Freshwater  Height,  and  these  values  were  taken  directly  from  the 
Depth  of  32 /oo  charts  of  freshwater  distribution.  The  valuer 

for  the  ratio  K  /K .  computed  in  this  manner  for  winter  observations  of 
Cruise  287  and  the  summer  observations  of  Cruise  290  ranged  from  10^  near¬ 
shore  to  10°  offshore.  Assuming  the  scale  of  the  phenomena  is  described 
by  the  r.m.s.  of  the  x  and  y  distances,  a  relationship  was  found  between  the 
eddy  coefficient  ratios  and  the  scale  of  the  phenomenon  (Figure  40: 

K/K  »  |2.6). 

y  - 
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If  K,  1  cm  sec  everywhere,  the  values  of  the  ratio  would  also  be 

values  for  "K  but  the  slope  K  :L  is  then  significantly  larger  than  has 
ever  been  reported.  If  this  is'^the  case  we  would  have  to  conclude  that  the 
simplified  equation  does  not  adequately  describe  the  effluent  distribution. 
However,  if  K  were  1  to  2  orders  of  magnitude  smaller  near  the  River  mouth 
than  offshore,  then  the  slope  of  the  line  in  Figure  40  would  approach  Kal8"''. 
This  latter  possibility  is  reasonable  in  view  of  the  general  decrease  in 
stability  of  the  water  column  with  increase  in  distance  from  the  mouth,  as 
K.  must  increase  with  a  decrease  in  stability.  In  either  event,  the  re¬ 
lationship  shown  in  Figure  40  is  so  consistent  that  it  can  be  used  in  an 
empirical  model  for  predicting  the  dispersion  of  the  plume  water  in  the 
offshore  environment. 

METHOD  II--The  horizontal  coefficients  were  evaluated  using  the 
Eulerian  expansion  of  the  distribution  of  variables  equation  assuming  a 
steady  state  and  no  vertical  diffusion.  The  assumptions  are  highly 
questionable  and  the  coefficients  calculated  are  too  large,  but  the  tech¬ 
nique  is  presented  to  show  the  order  of  magnitude  of  the  coefficients 
derived.  Data  from  Cruise  290  were  vised.  The  surface  salinities  from  that 
Cruise  (Fig.  19)  show  distinct  pockets  of  low  salinity  water.  An  assump¬ 
tion  was  made  that  these  cells  were  generated  from  the  river  in  a  similar 
manner  and  thus  were  originally  of  the  same  size.  As  they  move  away  from 
the  source  the  volume  increases  due  to  addition  of  surrounding  seawater. 

The  horizontal  gradient  of  salinity,  assumed  to  be  symmetrical  in  all 
directions,  was  determined  for  the  three  cells  located  approximately  along 
the  axis.  The  one  cell  remote  from  the  source,  but  close  to  the  coast, 
was  neglected  because  of  possible  reflection  from  the  coastal  barrier.  The 
rate  of  change  of  the  gradient  with  distance,  ,^2g  was  computed  by  find- 

dX? 

ing  the  difference  between  gradients  of  successive  cells  and  dividing  by 
the  distance  between  the  centers  of  the  cells.  Assuming  a  horizontal 
velocity  of  10  cm  sec" ^  as  a  reasonable  average  resultant  for  the  effect s 
of  the  geostrophic  and  jjind  driven  currents,  and  K.  were  calculated  to 
be  K  =  K  =  4  x  10 1  cm'  sec  and  the  age  of  the  most  remote  cell  to  be 
27  days.  ,vIf  the  single  low  salinity  cell  of  Cruise  2 88  (Fig.  18)  is  the 
same  as  the  most  remote  cell  in  Figure  19,  the  calculated  age  agrees  with 
the  actual  time  interval  between  Cruise  288  and  290  of  30  days. 


if  horiz 


From  the  above  two  approaches  it  appears  that  the  horizontal  dil'f  .  - 
sion  coefficients  are  of  the  order  of  10'  era*  sec”  •  with  a  l<  fii  It  1  tend 
in  the  offshore  regions  for  the  values  to  increase  with  distance  from  shor<  . 
In  the  nearshore  tidal  and  hydraulic  region  (less  than  0  km  from  the  rlv  ’• 
mouth)  the  coefficient  Is  considerably  Less.  This  area  is  now  being 
studied  in  more  detail. 

From  these  results  it  would  seem  that  an  empirical  model  could  be 
made  which  would  allow  prediction  of  the  distribution  of  the  Columbia  River 
plume.  By  p Lot  ting  the  weekly  trajectories  of  water  movement  it  is  possible 
to  arrive  at  a  predicted  axis.  Using  the  solution  of  Roberts  (1923)  and 
the  derived  eddy  coefficients  the  lateral  boundaries  of  the  plume  could  be 
predicted  from  month  to  month  throughout  the  year.  Application  of  this 
technique  to  the  dispersion  of  river  water  in  the  open  ocean  awaits  a  better 
evaluation  of  the  changes  of  the  eddy  coefficient  with  distance  from  the 
river  mouth.. 


SUMMARY 


The  freshwater  flow  from  the  mouth  of  the  Columbia  River,  located 
at  the  boundary  between  Oregon  and  Washington,  is  14  percent  of  the  annual 
runoff  from  the  United  States  and  represents  the  major  supply  of  freshwater 
to  the  open  ocean  along  the  coast  from  the  Strait  of  Juan  de  Fuca  to  San 
Francisco  Bay.  During  the  spring  the  peak  discharge  of  nearly  17,000  m-sec” 
(-.  (X  ,0C  c  ft  V  •  •"  ' )  is  associated  with  the  snow-melt  of  the  interior  mountain 
and  duritifc  winter  a  secondary  peak  of  about  7,000  m'sec"  1  (250,000  ft-Yec"^) 
corresponds  to  heavy  coastal  rainfall.  The  river  freshwater  has  a  signifi¬ 
cant  effect  on  biological,  geological,  chemical,  and  physical  character¬ 
istics  of  the  adjacent  open  sea.  After  the  river  water  is  mixed  witli  about 
2  parts  ambient  seawater  in  the  estuary  it  enters  the  nearshore  environment 
in  12-hour  pulses  with  each  ebb  tide  whe2-e  further  mixing  occurs  through 
hydraulic,  tidal,  and  surf  processes.  In  the  area  between  the  mouth  and 
70  to  50  km  seaward,  the  river  water  is  mixed  with  about  10  parts  of  ambient 
seawater  having  a  salinity  of  approximately  32.5$'-''.  Thence  the  effluent 
enters  into  the  offshore  ocean  regime  characterized  by  weak  north  setting 
currents  and  winds  in  the  winter  and  south  setting  currents  and  winds  in  the 
summer.  The  effluent  is  confined  to  the  surface  layer;  it  is  rarely  found 
below  4-0  meters  where  a  .harp  pycnocline  limits  downward  mixing  but  spreads 
Laterally  over  a  large  area.  In  the  ocean  regime  the  pancake- Like  effluent 
is  acted  upon  by  wind  stress,  geostrophic  currents,  and  turbulent  processes. 
During  the  stunmer  this  water  is  distributed  in  the  form  of  a  plume  extend¬ 
ing  southwest  almost  900  kilometers  from  the  river  mouth  to  near  40  N 
latitude.  During  the  winter  the  effluent  lies  in  a  40  to  55-km  wide  belt 
of  low  salinity  water  adjacent  to  the  coast  from  40  km  south  of  the  mouth  of 
the  Columbia  to  north  of  the  Strait  of  Juan  de  B\ica. 

This  description  of  the  dispersion  of  Columbia  River  water  in  the 
sea  is  based  on  12  oceanographic  surveys  conducted  during  UX'l  and  19''2 .  A 
technique  of  computing  wind  driven  transport  of  the  effluent  waters  has.  been 
developed  and  this  program,  along  with  detailed  observations  of  the  geos tr,  - 
phir  currents,  shows  promise  as  a  technique  for  predicting  the  axis  of  the 


plume.  Studies  of  the  lateral  and  vertical  diffusion  of  river  water  by 
turbulent  mixing  indicate  that  the  salt  flux  is  preferentially  upward  and 
the  vertical  eddy  coefficient  in  the  offshore  region  is  in  the  range  of  1 
to  10  cm^sec-!.  The  lateral  eddy  coefficient  was  estimated  by  two  tech¬ 
niques  and  found  to  be  about  lo”  cm^sec"^.  A  model  for  the  diffusion  of  the 
effluent  waters  in  the  offshore  region  could  be  developed  from  these  results 
but  application  of  the  model  awaits  the  evaluation  of  the  variation  of  the 
vertical  eddy  coefficient  with  distance  from  the  river  mouth. 


Future  Work 


It  is  realistic  to  divide  the  study  of  the  effects  of  the  Columbia 
River  water  in  the  sea  into  three  phases.  The  first  phase  is  the  study  of 
the  offshore  regions.  This  phase  was  the  major  concern  during  1961  and 
1962,  and  the  major  portion  of  this  study  is  reported  in  this  paper.  The 
second  phase  is  an  investigation  of  the  more  complicated  nearshore  environ¬ 
ment  from  the  river  mouth  to  the  limit  of  major  tidal  and  hydraulic  influ¬ 
ences.  Stations  spaced  closely  in  time  and  place  using  more  than  one  ship 
are  necessary  to  investigate  this  region  of  rapid  mixing.  Exploratory 
cruises  in  the  nearshore  regime  were  conducted  during  the  summer  of  1963  and 
plans  are  being  made  for  an  intensive  investigation  of  the  nearshore  region 
during  1964.  Emphasis  will  be  placed  on  mixing  rates,  circulation,  up- 
welling,  estimations  of  eddy  coefficients  and  a  general  three  dimensional 
description  by  multiship  sampling,  dye  studies,  moored  current  meter  studies 
and  radio  beacon  current  drogues. 

A  third  phase  involves  the  study  of  estuarine  circulation,  mixing, 
sedimentation,  and  biological  transport  of  suspended  and  dissolved  constitu¬ 
ents.  Past  investigations  have  been  limited  to  studies  by  the  U.S.  Army 
Engineers  in  1959>  oceanographic  exploratory  sampling  in  the  summer  of  i960 
by  the  University  of  Washington,  and  model  studies  by  the  U.S.  Army 
Engineers  and  U.S.  Geological  Survey. 

Estuary  studies  by  field  and  model  investigations  being  conducted 
by  others  will  lead  to  a  better  understanding  of  the  nearshore  and  offshore 
reg ime  s . 


The  investigations  envisioned  for  the  near  future  in  the  offshore 
area  are: 

1.  Studies  of  upwelling  and  divergence  areas. 

2.  Measurements  of  currents  near  the  Columbia  River  and 
Strait  of  Juan  de  Fuca. 

3.  Evaluation  of  eddy  coefficients. 

4.  Deep  and  bottom  current  measurements. 

5-  Detailed  studies  on  the  short-term  fluctuations  in 
mass  structure. 

6.  Observations  leading  toward  a  better  understanding 
of  baroclinic  response  times. 

7.  Total  transport  computations  through  an  extension  of 
the  wind  program. 
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the  above  addressee  to  the  British 
Joint  Services  Staff  for  further 
distribution  in  England  and  Canada.) 

1  Oceanographer 

Chief  ©f  Naval  Operations 
0P-09B6 

Washington  25,  D.  C- 
1  Commander 

U.  S j  Naval  Oceanographic  Office 
— ^Wasiyington  25,  D.  C. 

Library  (Code  1640 ) 

1  U.  S.  Naval  Branch 
Oceanographic  Office 
Navy  3923,  Box  77,  F.P.O. 

San  Francisco,  California 


Chief,  Bureau  of  Naval  Weapons 

Department  of  the  Navy 

Washington  25,  D.  C. 

1  Attn:  FASS 

1  Attn :  RU-222 


Nival 


Office  of  the  U.  S. 

Weather  Service 
U-  S.  Naval  Station 
Washington  25,  D.  C. 


Chief,  Bureau  of  Yards  &  Docks 
Office  of  Research 
Department  of  the  Navy 
Washington  25,  D.  C. 

At  tn :  Code  "JO 

Command ing  Officer  Sc  Director 
U.  S.  Navy  Electronics  Laboratory 
San  Diego  52,  California 
Attn:-  Code  2201 
Attn;  Code  2420 

Commariding  Officer  &  Director 
U.  S-  Naval  Civil  Engineering 
Laboratory 

Port  Huetieme,  California 
Attn-  Code  L54 

Code  31^5 
Box  7 

Pt.  Mugu  Missile  Range 
Pt.  Mugu,  California 

Commander,  Naval  Ordnance  Laboratory 
White  Oak,  Si Lver  Spring,  Maryland 
Attn:  E.  Liberman,  Librarian 


Command  mg  Off  ice 
Naval  Ordnance  rIV 
China  Lake,  Calif 
Attn:  Code  753 
Attn:  Code  508 


t  Station 


Chief,  Bureau  of  Ships 
Department  of  the  Navy 
Washington  25,  D.  C- 
1  Attn:  Code  373 

1  Of  l'lcer-  in-Cliarge 

U.  S.  Navy  Weather  Research 

Fac ility 

Naval  Air  Station,  Bldg.  R-l*8 
Norfolk,  Virginia 

1  U.  S  Fleet  Weather  Facility 
U.  S-  Naval  Station 
San  Diego  35,  California 

1  Commanding  Officer 

U.  S.  Navy  Air  Development 
Center 

Johnsville,  Pennsylvania 
Attn:  NADC  Library 

1  Supe  r intendent 

U.  S-  Naval  Academy 
Annapolis,  Maryland 

2  Department  of  Meteorology  & 

Oceanography 

U.  S.  Naval  Postgraduate  School 
Monte  ivy,  California 

1  Commanding  Officer 

U.  S-  Naval  Underwater  Sound 
Lnborat  ory 

New  London,  Connecticut 

1  Commanding  Officer 

U  S-  Navy  Mine  Defense  lab. 
Panama  City,  Florida 

1  Commanding  Officer 

U.  S-  Fleet  Weather  Central 
Navy  Department 
Washington  25,  D  C. 


1  Commanding  Officer  1  Off  ice r- in-Charge 

Naval  Rad-TT*+*igical  Defense  Laboratory  U-  S  Fleet,  Numerical  Weather 
San  Franc  isco7>0«NUf’orri  ia  Fae  ility 

Monterey,  California 

1  Commanding  Officer 

U.  S-  Naval  Underwater  Ordnance 
Station 

Newport,,  Rhixie  Island 


AIR  FORCE 


1  Hdqtrs.,  Air  Weather  Service 

(awss/tipd) 

U.  S.  Air  Force 

Scott  Air  Force  Base,  Illinois 

1  ARCRL  (CRZF) 

L.  G.  Hanscom  Field 
Bedford,  Massachusetts 


ARMY 

1  Array  Research  Office 

Office  of  the  Chief  of  R  &  D 
Department  of  the  Army 
Washington  25,  D.  C. 

1  U.  S.  Army  Beach  Erosion  Board 
5201  Li title  Falls  Road,  N.  W. 
Washington  16,  D.  C. 

1  Army  Research  Office 
Washington  25,  D.  C. 

Attn:  Environmental  Sciences  Div. 


OTHER  U.  S.  GOVERNMENT  AGENCIES 

1  Office  of  Technical  Services 
Department  of  Commerce 
Washington  25,  D.  C- 

20  Defense  Documentation  Center 
Cameron  Station 
Alexandria,  Virginia 

2  National  Research  Council 
2101  Constitution  Avenue 
Washington  25,  D.  C- 

Attn:  Committee  on  Undersea  Warfare 
Attn:  Committee  on  Oceanography 

1  Laboratory  Direst or 
Biological  Laboratory 
Bureau  of  Commercial  Fisheries 
P.  0.  Box  6121,  Pt.  Loma  Street 
San  Diego,  California 

1  Commandant  (OSR-2) 

U.  S.  Coast  Guard 
Washington  25,  D.  C- 


1  Commanding  Officer 

Coast  Guard  Oceanographic  Unit 
Bldg.  159,  Navy  Yard  Annex 
Washington  25,  D.  C. 

1  Director 

Coast  &  Geodetic  Survey 
U.  S.  Department  of  Commerce 
Washington  25,  D.  C. 

Attn:  Office  of  Oceanography 

1  Geological  Division 
Marine  Geology  Unit 
U.  S.  Geological  Survey 
Washington  25,  D.  C. 

1  Director  of  Meteorological  Research 
U.  S.  Weather  Bureau 
Washington  25,  D.  C. 

1  Director 

U.  S.  Army  Engineers  Waterways 
Experiment  Station 
Vicksburg,  Mississippi 
Attn:  Research  Center  Library 

1  laboratory  Director 

Bureau  of  Commercial  Fisheries 
Biological  Laboratory 
450 -B  Jordan  Hall 
Stanford,  California 

1  Bureau  of  Commercial  Fisheries 
U.  S.  Fish. &  Wildlife  Service 
Post  Off ic/e  Box  3830 
Honolulu  12,  Hawaii 

1  Attn:  Librarian 

1  Laboratory  Director 
Biological  laboratory 
Bureau  of  Commercial  Fisheries 
'P>Ns0.  Box  3098,  Fort  Crockett 
Galveston,  Texas 

1  Laboratory  Director 

Biological  Laboratory,  Auke  Bay 
Bureau  of  Commercial  Fisheries 
P.  0.  Box  1155 
Juneau,  Alaska 


1+ 


1  Laboratory  Director 
Biological  Laboratory 
Bureau  of  Commercial  Fisheries 
P.  0.  Box  6 

Woods  Hole,  Massachusetts 

1  laboratory  Director 
Biological  Laboratory 
Bureau  of  Commercial  Fisheries 
P.  0.  Box  280 
Brunswick.  Georgia 

1  Laboratory  Director 
Biological  Laboratory 
Bureau  of  Commercial  Fisheries 
P.  0.  Box  271 
La  Jolla.  California 

1  Bureau  of  Sport  Fisiieries  &  Wildlife 
U.  S.  Fish  and  Wildlife  Service 
Sandy  Hook  Marine  Laboratory 
P.  o'.  Box  428 
Highlands ,  New  Jersey 
Attn:  Librarian 

1  Director 

National  Oceanographic  Data  Center 
Washington  25.  D.  C. 

2  Defense  Research  Member 
Canadian  Joint  Staff 

2450  Massachusetts  Avenue,  N.  W. 
Washington  8,  D.  C 

2  Library,  U  S  Weather  Bureau 
Washington  25.  D  C- 

1  Director,  Biological  Laboratory 
Bureau  of  Commercial  Fisheries 
Navy  Yard  Annex 

Building  74 
Washington  25,  p.  C 

2  Director,  Bureau  of  Commercial 

Fisheries  \ 

U.  S.  Fish  &  Wildlife  Service 
Department  of  Interior 
Washington  25,  D  \C. 

1  Dr.  Orlo  E.  ChitStaA 

U.  S.  Geological  SunTey... 

345  Middle field  Road 
Menlo  Park.  California 


1  Dr.  John  S.  Schlee 

U.  S.  Geological  Survey 
c/o  Woods  Hole  Oceanographic 
Institution 

Woods  Hole,  Massachusetts 


RESEARCH  LABORATORIES 

2  Director 

Woods  Hole  Oceanographic  Inst. 
Woods  Hole,  Massachusetts 

3  Project  Officer 
Laboratory  of  Oceanography 
Woods  Hole,  Massachusetts 

1  Director 

Narragansett  Marine  Laboratory 
University  sf  Rhode  Island 
Kingston,  Rhode  Island 

1  Bingham  Oceanographic  Labs. 

Yale  University 

New  Haven,  Connecticut 

1  Gulf  Coast  Research  Laboratory 
Post  Office  Box 
Ocean  Springs.  Mississippi 
Attn:  Librarian 

1  Chairman,  Department  of 

Meteorology  and  Oceanography 
New  York  University 
New  York  53.  New  York 

1  Director 

Lamont.  Geological  Observatory 
Torrey  Cl iff 
Palisades.  New  York 

1  Director 

Hudson  Laboratories 
145  Palisade  Street 
Dobbs  Ferry,  New  York 

1  Great  Lakes  Research  Division 
Institute  of  Science  &  Technology 
University  of  Michigan 
Ann  Arbor,  Michigan 
I  Attn:  Dr.  John  C.  Ayers 


1  Dr.  Harold  Haskins 
Rut ge rs  Un ive rs  l ty 
New  Brunswick,  New  Jersey 

1  Director 

Chesapeake  Bay  Institute 
Johns  Hopkins  University 
121  Maryland  Hall 
Baltimore  18,  Maryland 

1  Mail  No.  J-3009 
The  Martin  Company 
Baltimore  3,  Maryland 
Attn:  J.  D.  Pierson 

1  Mr.  Henry  D.  Simmons,  Chief 
Estuaries  Section 
Waterways  Experiment  Station 
Corps  of  Engineers 
Vicksburg,  Mississippi 

1  Director,  Marine  Laboratory 
University  of  Miami 
#1  Rickenbacker  Causeway 
Virginia  Key 
Miami  1+9,  Florida 

1  Nestor  C.  L.  Granelli 
Department  of  Geology 
Columbia  University 
Palisades.  New  York 

2  Head,  Department  of  Oceanography 

&  Meteorology' 

TY'xas  A  &  M  CoiLege 
College  Stat ion,  Texas 

1  Director 

Scripps  Institution  of  Oceanography 
La  Jolla,  California 

1  Allan  Hancock  Foundat ion 
University  Park 
Los  Angeles  7,  California 

1  Head,  Department  of  Oceanography 
Oregon  State  University 
Corvallis,  Oregon 

1  Department  of  Engineering 
University  of  California 
Berkeley,  California 


1  Director 

Arctic  Research  Laboratory 
Barrow,  Alaska 

1  Dr.  C-  I.  Beard 

Boeing  Scientific  Research  Labs. 

P.  0.  Box  3981 
Seattle  2k,  Washington 

1  Head,  Department  of  Oceanography 
University  of  Washington 
Seattle  5,  Washington 

1  Geophysical  Institute  of  the 
University  of  Alaska 
College,  Alaska 

1  Director 

Bermuda  Biological  Station 
for  Research 
St.  Georges,  Bermuda 

1  Technical  Information  Center,  CU-201 
Lockheed  Missile  and  Space  Division 
3251  Hanover  Street 
Palo  Alto,  California 

1  University  of  Pittsburgh 
Environmental  Sanitation 
Department  of  Public  Health  Practice 
Graduate  School  of  Public  Health 
Pittsburgh  13,  Pennsylvania 

1  Director 

Hawaiian  Marine  Laboratory 
University  of  Hawaii 
Honolulu,  Hawaii 


Dr.  F.  B.  Berger 
General  Precision  Laboratory 
Pleasantv ille ,  New  York 

Dr.  J.  A.  Gast 
\  Wildlife  Building 

Humboldt  State  College 
Areata,  California 

lj  Department  of  Geodesy  &  Geophysics 
|  Cambridge  University 
j  Cambridge,  England 


1  Applied  Physics  Laboratory 
University  of  Washington 
1013  NE  Fortieth  Street 
Seattle  5,  Washington 

1  Documents  Division  -  ml 

University  of  Illinois  Library 
Urbana,  Illinois 

1  Director 

Ocean  Research  Institute 
University  of  Tokyo 
Tokyo,  Japan 

1  Marine  Biological  Association 
of  the  United  Kingdom 
The  laboratory 
Citadel  Hill 
Plymouth,  England 

1  ASW  Information  Research  Unit 
Building  80,  Plant  A-l 
Lockheed-California  Company 
Burbank,  California 


1  Institute  of  Geophysics 
University  of  Hawaii 
Honolulu  l4,  Hawaii 

1  Dr.  Wilbur  Marks 
Oceanics,  Inc. 

114  East  40th  Street 
New  York  16,  New  York 

1  Mr.  Neil  L.  Brown 

Bissett -Berman  Corporation 

"G"  Street  Pier 

San  Diego,  California 


1  New  Zealand  Oceanographic  Institute 
Department  of  Scientific  and 
Industrial  Research 
P.  0.  Box  8009 
Wellington,  New  Zealand 
Attn:  Librarian 


1  President 

Osservatorio  Geofisico  Sperimentale 
Trieste,  Italy 


1  Advanced  Research  Projects  Agency 
Attn:  Nuclear  Test  Detection  Office 
The  Pentagon 
Washington  25,  D.  C. 


1 


